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Abstract. A special class of integrable nonlinear differential equations related to A. Ill-type 
symmetric spaces and having additional reductions are analyzed via the inverse scattering 
method (ISM). Using the dressing method we construct two classes of soliton solutions 
associated with the Lax operator. Next, by using the Wronskian relations, the mapping 
between the potential and the minimal sets of scattering data is constructed. Furthermore, 
completeness relations for the 'squared solutions' (generalized exponentials) are derived. 
Next, expansions of the potential and its variation are obtained. This demonstrates that 
the interpretation of the inverse scattering method as a generalized Fourier transform holds 
true. Finally, the Hamiltonian structures of these generalized multi-component Heisenberg 
ferromagnetic (MHF) type integrable models on A. Ill-type symmetric spaces are briefly 
analyzed. 
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1 Introduction 

It is well known [9], that two important classes of integrable nonlinear evolution equations 
(NLEE) are related to symmetric spaces: the multi-component nonlinear Schrodinger equations 
and their gauge equivalent Heisenberg ferromagnet (HF) equations: 

iSt = ^[S,Sxx], S'^{x,t) = l. 

Here S takes values in some simple Lie algebra. In the simplest case of the algebra su(2) the 
function S describes the spin of one-dimensional ferromagnet [4]. 

The HF's equation is integrable in the sense of inverse scattering transform [9, 45]. Its Lax 
representation has the form: 

L(A) = id., - XS, A{X) = idt + '-^ [S, 5J + 2X^S. 

Since the time the complete integrability of HF equations was discovered, many attempts for 
constructing its generalization have been made [22]. The main goal of this paper is to derive 
and analyze a special class of NLEE which is obtained from the HF type equations by applying 
additional algebraic symmetries. 

A well known method for constructing new integrable NLEE is based on the so-called re- 
duction group, introduced in [32, 33, 34, 35] and further developed in [13, 30, 28, 29]. It led to 
the discovery of the 2-dimensional Toda field theories [34, 36]. Its latest developments of the 
method led to the discovery of new automorphic Lie algebras and their classification [30, 28, 29]. 

The hierarchies of integrable NLEE are generated by the so-called recursion (generating) 
operator. Such operators have been constructed and analyzed for a wide class of Lax operators L 
and appeared to generate not only the Lax representations, but also the hierarchy of NLEE's 
related to a given Lax operator L, their conservation laws and their hierarchy of Hamiltonian 
structures, see [7, 40, 45, 20, 21] and the numerous references therein. Such operators can 
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be viewed also as Lax L operators, taken in the adjoint representation of the underlying Lie 
algebra 3. 

The (generating) recursion operator appeared first in the AKNS-approach [1] as a tool to 
generate the class of all ^d-operators as well as the NLEE related to the given Lax operator L{\). 
Next, LM. Gel'fand and L.A. Dickey [10] discovered that the class of these A-operators is 
contained in the diagonal of the resolvent of L. The kernel of the resolvent of L can be explicitly 
constructed in terms of the so-called fundamental analytic solutions (FAS), see [11, 17, 12]. 

The construction of the recursion operator for Lax operators, whose explicit dependence on 
the spectral parameter A is comparatively simple (say, linear, or quadratic) was done a long time 
ago [1, 27, 15, 16]. Furthermore, the completeness property for the set of eigenfunctions of the 
recursion operator (the 'squared solutions' of L) is of paramount importance. The complete- 
ness of the 'squared solutions' plays a fundamental role in proving that the inverse scattering 
method is, in fact, a nonlinear analogue of the Fourier transform, which allows one to linearize 
the NLEE. Using these relations one is able to derive all fundamental properties of the NLEE 
on a common basis. 

In the symmetry approach to integrable equations the so-called 'formal recursion operator' 
plays a crucial role. It has the general property to map a symmetry into a symmetry of the given 
integrable NLEE [41]. If the considered NLEE possesses an infinite hierarchy of symmetries [26], 
or conservation laws [43] of arbitrary (high) order, or can be linearized by a differential substi- 
tution [44], then it has a formal recursion operator [39, 38, 42]. For an extensive and up-to-date 
surveys on symmetry approach to integrability we refer to [2, 37] and the references therein. 
We just note that in the Symmetry approach the formal recursion operator does not depend on 
the boundary conditions, imposed on the potential of the associated Lax operator. An alterna- 
tive approach for construction of recursion operator from a given Lax operator is presented by 
M. Giirses, A. Karasu and V. Sokolov in [24]. This method is pure algebraic, it is not related 
yet with the spectral theory of the Lax operator L and it is not sensitive to the choice of the 
class of admissible potentials of L. 

The problem of deriving recursion operators becomes more difficult when we impose ad- 
ditional reductions on L. If this additional reduction is compatible with L, being linear or 
quadratic in A, the construction of the recursion operators is not a difficult task (see [20]; an 
alternative construction of A is given in [24, 23]). The effect of the Z^-reduction is as follows: 
i) the relevant 'squared solutions' have analyticity properties in sectors of the complex A-plane 
closing angles vr/n; ii) the grading of the Lie algebra g = (B^Z^Q^^^^ is more involved and as 

n-l 

a consequence the recursion operator is factorized into a product of n factors A = A^, and 

k=Q 

each of the factors A^ maps A^: g^^"^) — )■ 

The applications of the differential geometric and Lie algebraic methods to soliton type equa- 
tions led to the discovery of a close relationship between the multi-component (matrix) inte- 
grable equations (of nonlinear Schrodinger type) and the symmetric and homogeneous spaces [9] . 
Later on, this approach was extended to other types of multi-component integrable models, 
like the derivative NLS, Korteweg-de Vries and modified Korteweg-de Vries, N-wave, Davey- 
Stewartson, Kadomtsev-Petviashvili equations [3, 8]. 

The purpose of the present paper is to derive nonlinear evolution equations to generalize 
Heisenberg's model and study some of the properties of their Lax operators. We are going to 
focus our attention on equations whose Lax representation is related to su(3). 

This paper is a natural continuation of our previous papers [18, 19, 14]. In Section 2 below 
we give some of the necessary preliminaries. We also study the Z2-reductions of the generalized 
Heisenberg ferromagnets (Z2-HF) related to symmetric spaces and outline the spectral pro- 
perties of the relevant Lax operator and the construction of its fundamental analytic solutions. 
Section 3 is dedicated to the soliton solutions of the (Z2-HF) models. In their derivation we make 
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use of the dressing method [40, 49, 50, 51]. Due to the additional reductions the Lax operator L 
may possess two types of discrete eigenvalues: generic ones, coming in quadruplets ±Afc, ±A^ 
and purely imaginary ones coming as doublets ziziKj. Therefore we will have two types of 
soliton solutions which we will call quadruplet and doublet solitons respectively. We outline 
the purely algebraic construction for deriving the A'^-soliton solutions for both types of solitons 
and provide the explicit expressions for the one-soliton solutions. Section 4 is devoted to the 
derivation of the recursion operators A. Here we first derive A using the Giirses-Karasu- 
Sokolov (GKS) method [24]. Next we analyze the Wronskian relations as a basic tool in the 
inverse scattering method [5, 6]. From them, there naturally arise the 'squared solutions', which 
play a fundamental role also in the analysis of the mapping between the set of admissible 
potentials and the minimal sets of scattering data. The 'squared solutions' can also be viewed 
as eigenfunctions of the recursion operators A± , which can be used as an alternative definition 
of the recursion operators. One can check that both approaches lead to equivalent expressions 
for A-I-. Section 5 is dedicated to the spectral properties of the recursion operators. There we 
start by proving the completeness relation for the 'squared solutions'. Next we use this relation 
to derive the expansions of ad^^^Li^^^, ad^^^ L2,x and adj^^ 6Li over the 'squared solutions'. In 
the last Section 6 these expansions are treated as generalized Fourier transforms, allowing one to 
linearize the NLEE. We also demonstrate that all fundamental properties of the class of NLEE 
can be formulated through the recursion operators. These include not only the description of 
the class of NLEE related to L, but also their integrals of motion and hierarchy of Hamiltonian 
structures. We end by discussion and conclusions. In Appendix A we have collected some useful 
formulae for the Gauss factors of the scattering matrix and explicitly derive the operator ad^^. 



2 Preliminaries 

2.1 Basic notations and general form of equations 

The main object in our paper is the following 2-component system of NLEEs: 

iut + Uxx + {uu*x + vvl)ux + {uul + vvDxU = 0, 

ivt + Vxx + {uu*x + vvl)vx + {uul + vvDxV = 0, (2.1) 

where the functions u : M x M — )• C and u : M x M — )■ C are assumed to be infinitely smooth and 
satisfy the following boundary conditions 

lim u{x,t) = 0, lim v{x,t) = e''^^ , (/)± G M. (2.2) 

X— ^-itoo X— >itoo 

Moreover, u and v are not functionally independent, but obey the constraint |np + |f p = 1. 

The system (2.1) represents a reduction of generalized Heisenberg ferromagnet equations 
related to the symmetric space SU{3)/ S{U{1) x U{2)), see [18, 19]. It possesses a zero curvature 
representation with Lax operators in the form 

L{X) =idx + XLi{x,t), (2.3) 
A{X) = idt + A^i(x, t) + X^A2ix, t), (2.4) 

where A is a spectral parameter and the matrix coefficients (potentials) read 

Li = I n* I , A2 = ^IL-L? = -| - 2/3 u*v |, (2.5) 
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^1 = a* , a = iux + i{uu*^ + vv*)u, b = ivx + ^(uu* + vv*)v. (2.6) 



The specific structure of the matrices above is a result of the simultaneous action of two Z2 
reductions on generic Lax operators L and A, namely 

Lt(A*) = -L(A), A^X*) = -A{X), (2.7) 
CL(-A)C = L(A), CA{-X)C = A{\), C = diag(l, -1, -1), (2.8) 

where the operation" is defined as follows 

L{X)ip{x, t, A) = idxip{x, t, A) — Xip{x, t, X)Li{x, t, A). 

The latter reduction represents Cartan's involutive automorphism [25, 31] involved in the defi- 
nition of the symmetric space SU{3)/ S{U{1) x U{2)), that is, it induces a Z2-grading in the Lie 
algebra s[(3, C) 

s[(3) = 5[°(3) © sl^(3), 5r (3) = {X Gg\ CXC = {-ITX}. 

It is evident that the grading condition 

[sr(3),sr(3)] Csr+^('^°'^2)(3) (2.9) 

is fulfilled. This grading will be used widely in our further considerations. Due to the existence 
of the grading (2.9), any function X{x, t, A) with values in s[(3) can be represented in the form: 

X{x,t,X)=X^{x,t,X)+X\x,t,X), X°'i(x,t,A) G s[°'\3). (2.10) 

Each component X^'^, in its turn, splits into a term commuting with Li and its orthogonal 
complement, namely 

= X^'^ + L2 = -A2 = Ll-\, (X0'^,L2)=0, (2.11) 

X^ = X^'-^ + KiLi, (X^'-^, Li) = 0. (2.12) 

The brackets above stand for Killing form, defined as: 

{X,Y) =ii{XY). 

It is evident that the functions kq and ki can be expressed as follows 



(X0,L2) {X\L 



{L2,L2)' ' {Lull}- 
One can easily check that the norms of Li and L2 are 

(Li,Li)=2, (L2,L2) = ^ 

and hence kq and ki are given by 

Ko = ^(XO,L2), Ki = ^{X\Li). 
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As we shall convince ourselves in some cases it proves to be more convenient to deal with Lax 
operators 

L{X)=id^ + Uo{x,t) + XJ, J = diag(l,0,-l), (2.13) 
A{X) = idt + Ao{x,t) + XAi{x,t)-X'^I, / = diag(l/3, -2/3, 1/3) (2.14) 

gauge equivalent to (2.3), (2.4). The gauge transformation, which puts Li and A2 into their 
diagonal form, is given by the unitary matrix 




Further on we are also going to use the explicit form of Uq given by next formula 



UU* + VV* \pi{uVx — VUx) uu* + vv* 



Uo = i9-^gx = -{ -V2iu*v*-v*ul) 2{u*ux + v*vx) -V2{u*vl - v*ul) \. (2.15) 
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«n* + vv* \f2{uvx — VUx) uu* + vv^ 



2.2 Direct scattering problem for the Lax operator (2.3) 

Here we shall provide a short summary of some basic results and notions concerning the spec- 
tral theory of the Lax operator L, the direct scattering problem and introduce the so-called 
fundamental analytic solutions. The spectral properties of L depend on the choice of the class 
of admissible potentials, i.e. the potentials are a subject to certain boundary condition. In the 
case of boundary conditions of type (2.2), the continuous spectrum of L coincides with the real 
axis in the complex A-plane, see [19]. 

In order to formulate the direct scattering problem for L, one needs to introduce fundamental 
sets of solutions^ Tp to the auxiliary (spectral) linear system 

L{X)'il){x,t,X) = idxtp{x,t,X) + XLi{x,t)'il){x,t,X) = 0. (2.16) 
Since L{X) and A{X) commute, fundamental solutions ip also satisfy the equation 

A{X)'ip{x, t, A) = idt^ix, t, A) + {XAi{x, t) + X^A2ix, t))'ip{x, t, A) = V(x, t, A)/(A). (2.17) 
The matrix-valued function 

/(A)= lim g^\XAi{x,t) + X^A2{x,t))g± = -X^I (2.18) 

I'— >±oo 

is called dispersion law of the nonlinear equation. 

A special type of fundamental solutions are the so-called Jost solutions -i/'i which are nor- 
malized as follows 

-iXJx „— 1 



hm V±(x,t,A)e-^-^-^^5i =1, 

X— >±oo 



where 



1 -1 

g±= lim g{x,t) = -= V2e'^± 



^F\irther we shall simply call them fundamental solutions for short. 
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diagonalizes the asymptotics Li ± = lim Li{x, t). Due to (2.18) one can show that the asymp- 

' x—>-±oo 

totic behavior of ^^J± do not depend on t, i.e. the definition is correct. The transition matrix 

r(t,A) = [iP+{x,t,X)]-'^ij.{x,t,X). (2.19) 

is called scattering matrix. It can be easily deduced from relation (2.17) that the scattering 
matrix evolves with time according to the linear differential equation 

idtT+[f{\),T] = 0, 

which is integrated straight away to give 

T{t,X) = e*-^W*T(0,A)e-^-^(^)*. 

Since the time parameter t does not play a significant role in our further considerations, we shall 
omit it (it will be fixed). 

The action of Z2-reductions (2.7), (2.8) imposes the following restrictions 

[i^i{x,X*)r' = ^±{x,X), [rt(A*)]-' = r(A), (2.20) 
CV'±(x,-A)C = ^±(a;,A), CT(-A)C = r(A) (2.21) 

on the Jost solutions and the scattering matrix. 

From now on we will assume that ^+ = 0. Hence the asymptotic values of Li and g 

become equal to each other: 

Li,± = Ll,as, 9+ = 9- = 5as- 

The complex A-plane is separated by the continuous spectrum of L (real axis) into two regions 
of analyticity: the upper half plane C+ and the lower one C_ . We are going to sketch two ways 
of constructing fundamental solutions (FAS) x~^{x, A) and x~{xj ^) which are analytic functions 
in C4. and C_ respectively, see [19] for more detailed explanations. The first way is based on 
introducing some auxiliary functions 

77±(x,A)=5a-sV±(a;,A)e-^^■^^ 
which satisfy the auxiliary system: 

idx'n± + U{x,X)rj±{x,X) - Xr]±{x,X)J = 0, U{x,X) = Xg^^Li{x)gas 

with the boundary conditions lim r]±{x,X) = U. Equivalently, r]±{x,X) can be regarded as 

a;— > it 00 

solutions to the following Volterra-type integral equations: 



77±(x,A) = U + i r dye^^-^("-s^)[C/(y,A)-AJ]77±(y,A)e-^^^("-s^). 

J±oo 

Then we define ^~^{x, A) as a solution to the following set of integral equations: 

CU^,X) = SM + i r (iye*^^-'"-^")^"-^) [(f/(y,A)-AJ)e+(y,A)],^, k<l, 
J —00 

?+(x. A) = i rdye^A(j,,(A)-j„)(.=-i/) [(c/(y, A) - AJ)e+(y, A)],, , k>l, 
J 00 

and ^~(x, A) being a solution to 

4-;(x, A) = dki+i r dye*^(^'='=-^")(--?/) {(U{y, A) - XJ)r{y, A)),, , k<l, 
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X 




Q(a:,A) = i / dye*^(-^--^")(^-J')((C/(y,A)-AJ)r(y,A)),,, k>l. 

J — oo 

It is easy to check that ^+ and have the analytic properties in C+ and C_ respectively 
due to the appropriate choice of the lower integration limits. The initial fundamental analytic 
solutions x^(3^)A) are obtained from ,^^(x,A) by applying the transformation: 

X^{x,X)=g.st{x,X)e'''^. (2.22) 

Another way of how one can construct FAS is by multiplying the Jost solutions by appropriately 
chosen matrices. It turns out that these matrix factors are involved in the Gauss decomposition 

r(A) = T^{X)D^{X){S^{X))-^ (2.23) 

of the scattering matrix T(A). Here is a list of all Gauss factors in a parametrization we are 
going to use further in our exposition: 

S^{X)= exp (A) , r± (A) = exp r± (A) , 

.+ (A) = I I , r+(A) 





s-{X)= { s^^ I , r-(A) 

D'''(A) = diag(mj['', m^/rri]^, 1/m^), D~(A) = diag(l/m2 , m2 m^), (2.24) 

where m^(A) (resp. m^(A)) are the principal minors of T(A) of order k = 1,2. Then x'^ and x~ 
are expressed as follows 

X^ix, A) = ^_(x, X)S^ = V'+(x, X)T^iX)D^{X). 

The relation above can be rewritten in the following manner, using (2.22): 

^+{x,X) = r{x,X)G{x,X), G{x,X) = e'^-^''{S-)-^S+{X)e-'^-^^. (2.25) 

Thus FAS can be regarded as solutions to a local Riemann-Hilbert problem. The established 
interrelation between the inverse scattering method and local Riemann-Hilbert problem proves 
to be useful in constructing solutions to NLEEs. 

It can be shown that the reduction conditions (2.20), (2.21) and equation (2.23) lead to the 
following demands on the Gauss factors 

(s+(A*))t = {s-{x)y\ cs^{-x)c = s^{x), 

(T+(A*))t = (r-(A))-\ CT±(-A)C = T^(A), C = g^,Cg-,' 

(i?+(A*))t = {D-{X))-\ CD^{-X)C = D^{X). 

According to (2.24) these relations can be also written as: 

(s+(A*))t = -.-(A), (r+(A*))t = -r-(A), 
Cs+(-A)C = s-(A), Cr+(-A)C = r-(A). 

Finally, combining all this information we see that the FAS obey the symmetry conditions 

ix^)\x,X*) = [x~ix,X)]-\ Cx+(x,-A)C = x~ix,X). 
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Remark 1. The Riemann-Hilbert problem allows singular solutions as well. The simplest types 
of singularities are simple poles and zeroes of the FAS and generically correspond to discrete 
eigenvalues of the Lax operator L. Due to the reduction symmetries the discrete eigenvalues 
must form orbits of the reduction group. Generic orbits contain quadruplets, so if /i is an 
eigenvalue, then and ±/i*, are eigenvalues too. However, we can have degenerate orbits too. 
If the eigenvalue lies on the imaginary axis we will have doublets of eigenvalues. 



3 Dressing method and soliton solutions 

In the present section we are going to derive and analyze the 1-soliton solution to the 2- 
component system (2.1). For this to be done, we are going to apply the dressing method 
proposed in [49] and developed in [50, 51, 34, 35]. In the previous section (see Remark 1) we 
established that the operator L may possess two types of discrete eigenvalues: one, coming in 
quadruplets as well as purely imaginary ones coming as doublets ziziKj. Therefore we will have 
two types of soliton solutions: generic (or quadruplet) ones and doublet ones (kinks). 

3.1 Rational dressing 

The idea of the dressing method is the following. Firstly, one starts from known solutions 
uo{x,t), vo{x,t) of (2.1) and a solution 'ipo{x,t,X) of the auxiliary linear problems 

lW(a)V'o = ^^.■V'o + a4°Vo = 0, 

A(o)(A)Vo = idt^o + {XA^i^ + X'A^2^)i^o = 0, (3.1) 

where L^i \ A^'^ and A'^^ have the same form as given in (2.5), (2.6) but u and v are replaced by 
uq and vo respectively. Then one constructs another function 'tpi{x,t,X) = ^{x,t, X)^o{x,t, X) 
which is assumed to be a fundamental solution to a similar system of equations 

L«(A)V'i = i^.V'i + a4'Vi = 0, 

yl(i)(A)V'i = idti'i + (A^S^) + A24'Vi = 0- (3-2) 

The matrices A^^^ and A^'^ have the same structure as Lf'\ Af'^ and A^^ but instead of 
-^0(3^,^) and vo{x,t) one plugs some new solutions ui{x,t) and vi{x,t) to (2.1) which are to be 
found. 

The dressing factor <I>(x,t,A) is assumed to be regular at |A| — >• 0, 00. Due to the reduction 
conditions (2.7), (2.8) the dressing factor has the symmetries: 

C<l>{x,t,-X)C = <^{x,t,X), (3.3) 
^{x,t,X)<^\x,t,X*) = t. (3.4) 

From (3.1) and (3.2) it follows that ^{x,t,X) is a solution to the following equations: 



ida:^ + Xl'ih - X<^Lf^ = 0, (3.5) 
idt^ + (AA;^) + A^A^^))^ - $(AAf) - A^Af ) = 0. (3.6) 

Since $(x, t, A) is regular at | A| -^00 from (3.5) one can derive the following relation between L^^^ 
0. 

L'^^\x,t) = ^{x,t,oo)Lf\x,t)^\x,t,oo). (3.7) 



and 4°^: 
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This equation will play a fundamental role in our further considerations since it allows one to 
generate a new solution to (2.1) from a given one. 

If we consider a A-independent dressing factor then from (3.5) and (3.6) we can deduce that $ 
do not depend on x and t either. Thus we have a simple rotation, related to the U{2) symmetry of 
the model. In order to get something more nontrivial we shall assume that ^{x, t, A) is a rational 
function of A with a minimal number of simple poles. From the reduction condition (3.3) it 
follows that if /i is a pole of $ then — /U is a pole too. It is natural to assume that these poles do 
not overlap with the continuous spectrum of L, i.e. Im n ^ 0. On the other hand (3.4) leads to 
the conclusion that has poles at /i* and — /i*. 

We first consider the generic case when fi ^ — /U*, so the poles of $ do not coincide with the 
poles of <I>~^ and choose the following anzatz for the dressing factor and its inverse: 

^ XM(x,t) \CM(x,t)C 

^{x,t,X) = TL+ , ' ^ + , , , (3.8) 

A — /X A + ^ 

^ w XX ^ XMUx,t) XCMUx,t)C 
X — fj,* X + fi* 

From (3.8) it follows that the matrix M satisfies the following equation: 

If we assume that the residue M is an invertible matrix (rank M = 3), then it is proportional to 1 
and the dressing is trivial. In order to obtain a non-trivial dressing M ought to be singular. For 
our purposes it suffices to consider the case rankM = 1. Then the residue can be represented 
as a product M = |n)(m| of a vector |n) = (ni,n2,n3)"^ and a co-vector (m| = {ml, 1712, m'^). 
Substitution of this representation into (3.9) leads to a linear equation for the vector |n): 

^ /x1n)(m|m) ^ /^1Cn)(m|Cm) ^ ^ 

Introduce the functions 

/o(2;,t) = -(m(x,t)|m(x,t)), /i(x,t) = -(m(x,t)|C|m(x,t)), ^ = uj + iK. 

Then, the solution of (3.10) reads: 

I f _ 2(i/i(x,t)|m(x,t))-/o(x,0|C|m(x,t))) 

The vector |m) is an element of the projective space CP^. Indeed, it is easily seen from above 
that a rescaling |m) — )• /|m) with any nonzero complex / does not affect the matrix M. 
Taking into account the anzatz (3.8) one can rewrite (3.7) as: 

= (1 + M + CMC)Lf\t + M + CMC)t. (3.12) 

Notice that the dressing preserves the matrix structure of the Lax operator L, since the dressing 
factor U + M + CMC is a block-diagonal matrix. From this equality and from (3.12) it follows 
that the potentials ui, vi can be expressed by no, fo in the following way: 
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where 

G = a;(m|m) + iK(m|C|m), = |/ip|mi|^ + (/i*)^|m2p + /U^lmsp, F = AiuJKm^m2. 

It is easy to verify that the matrix M is unitary for a generic choice of /x and the nonzero 
vector |m). 

Thus we expressed all quantities needed in terms of |m). It remains to find |m) itself. For 
that purpose we rewrite equations (3.5), (3.6) in the form: 

$(x,t,A)(i5^. + ALf^)$-^(x,t,A) = ALf\ 

^x,t,\){idt + XAf^ + X^A'^^^)^~^{x,t,\) = XA^l^ + \^A^^\ (3.14) 

It is obviously satisfied at A = 0. From (3.14), it follows that the residues at A = ±/x*, A = zt^u 
should vanish. It is sufficient to vanish the residue at A = ^u* (vanishing of the other residues 
follows from the symmetry conditions (3.3)). Computing the residue at A = //* and taking into 
account equation (3.9), we get 

V H- H- H' ' H' / 

i.e. |m(x,t)) is an eigenfunction of the dressed Lax operator. A general solution of equa- 
tions (3.15), (3.16) is 

|m(x,t)) = f{x,t)'ipo{x,t;fi*)\mo), 

where f{x,t) is an arbitrary non-vanishing complex function and |mo) G is a non-zero, but 
otherwise arbitrary complex vector. Without any loss of generality we can set f{x, t) = 1 (see 
the discussion above about the projective nature of the vector |m). 
Thus we proved the following proposition: 

Proposition 1. Let uo(x,t), VQ{x,t) form a solution of the system (2.1) and ^(){x,t,X) be 
a simultaneous fundamental solution of (3.1). Let also Re// ^ 0, Im/i ^ and |mo) G 
be a complex vector. Then ui(x,t),vi{x,t) given by (3.13) for mt being components of the 
vector |m) = 7/^0(2;, //*)|mo) satisfy (2.1) as well. The corresponding solution ^i(x,f. A) of the 
associated linear system (3.2) is given by tpi = ^{x,t, X)ipo{x,t, X) where $(x,t,A) is defined 
by (3.8), (3.11), (3.15) and (3.16). 

The new solution ui(x,t), vi{x,t) of equations (2.1) and the fundamental solution ipi(x,t; A) 
of the corresponding linear system are parameterized by a complex number /i and a complex 
vector I mo) G C'^. 

Let us now consider the special case when fj, = in. Then the sequence of steps necessary to 
determine $ slightly changes. Indeed, suppose we have a dressing factor in the form: 

*(..*,A) = 1 + a(QM + 5™). (3^,V) 

\ \ — in \ + lK J 

Due to the reduction (3.4) its inverse has the same poles as ^ itself and therefore the equation 
= U already contains second order poles. Vanishing of the poles of second and first order 
respectively leads to the following algebraic conditions: 



PCP^ = 0, 
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(^1 + P+ CP+C + p(^t + CP^C + = 0. (3.18) 

P is a degenerate matrix which means that it can be presented as P = |q)(p|. Then the former 
algebraic constraint transforms into a quadratic relation for the vector |p) 

(p|C|p) = 0. (3.19) 

Equation (3.18) is reduced to a linear system for |q) 

i + Sf^)c|p,...|<i. 

where a is some auxiliary real function to be found. The above linear system allows one to 
express |q) through |p) and a 

ia-^^^C^ C|p). (3.20) 

In order to find |p) and a we consider again the partial differential equations (3.14). The require- 
ment that the poles of second and first order in (3.14) vanish identically yields to differential 
relations for |p) and a: 

id^{p\ - iu{p\L^^^ = 0, idt{p\ - (p| (iz/4°^ + M^Af^) = 0, 

id^a = K(p|Lf^C|p), idta = k(p| (4°^ + 2ii/4°^)C|p). 

It is not hard to verify that \p{x,t)) and a{x,t) are expressed through the initial solution 
^po{x,t,X) as follows: 

\p{x,t)) = ipo{x,t,-iK)\po), (3.21) 

a{x, t) = -k{po\tP^^{x, t, iK)7po{x, t, m)C|po) + ctq, (3.22) 

where |po) and cjo are constants of integration and tJjq := 9a "00- 

After substituting (3.21) into (3.19) and taking into account the first Z2 reduction we see 
that the components of the polarization vector |po) are no longer independent but satisfy the 
constraint: 

(po|C|po>=0 ^ KiP = |po,2p + |po,3|^. (3.23) 

Thus to calculate the soliton solution itself one just substitutes the result for |p) and (p| 
into M and uses formula (3.12). Now we are able to formulate a result quite similar to Propo- 
sition 1: 

Proposition 2. Let it be given functions uo{x,t), vo{x,t) to satisfy (2.1) and 'ijjQ{x,t, X) to be 
a fundamental solution of (3.1). Let also k S ]R\{0}, cto G M and |po) he a non-vanishing 
complex vector satisfying (3.23). Then ui{x,t), vi{x,t) given by 

ui \ ^ IffiP + / icr + |p2p - |P3P 2p^P3 ^ f "° 

^1 J {\Pi\^ - \ '^P2P*3 icr - \P2\^ + \P3\^ J \ vq 

wherepk are the components of the vector |p) = ■iljQ{x,t, — zk)|po) is a solution of the system (2.1) 
too. The solution ^/;i(x,t, A) of the linear system (3.2) is given by -01 = ^ipo where <I> is defined 
by (3.17), (3.20), (3.21) and (3.22). 

As it is seen the new solution is parametrized by the polarization vector |mo), the real 
number ctq and the pole in. 

One can apply the dressing procedure repeatedly to build a sequence of exact solutions 
to the system which is parametrized by a set of complex constants fi^ and constant vectors 
|mfc) S CP^ (resp. a set of real constants kj, aoj and complex three- vectors |pfc) in the case of 
purely imaginary poles). 
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3.2 One-soliton solutions 



Let us apply the dressing procedure to a trivial solution uq = 0, vq = 1 of equation (2.1). In 
this case 



ipo{x,t,X) 



^ cos(Ax) exp 


^ i sin(Aa;) exp 



3 



3 



exp 





3 





i sin(Ax) exp 


cos(A3;) exp 



iX?t 
3 



(3.24) 



We are going to consider the generic case first, i.e. we have 4 distinct poles of $ and ^>~^ to form 
a 'quadruplet' {^u, — /i*}. It is convenient to decompose a constant complex vector |mo) 

according to the eigen-spaces of the endomorphism ■0o (3.24): 



|mo; = a |oj+/3foj+7fl 



(3.25) 



where q, /3, 7 are arbitrary complex constants. 

If the vector |mo) is proportional to one of the eigenvectors of the endomorphism V'O) then 
the corresponding matrix M does not depend on the variables x and t (due to the projective 
nature of the vector |m)) and the corresponding solution (3.13) is a simple unitary rotation of 
the constant solution tio = 0, = 1. 

Elementary solitons correspond to vectors |mo), belonging to an essentially two-dimensional 
invariant subspaces of the endomorphism -00 • In other words, elementary soliton solutions 
correspond to the choices of vector |mo) with only one zero coefficient in the expansion (3.25). 
Let us consider each of these three cases in more detail. 



Case (i): 07^0, /Jt^O, 7 = 

In this case, the solution does not depend on the variable t. It follows from (3.13) that the 
solution can be written in the following form 

Ml = 0, (3.26) 

_ / u;(|«pe^^^ + |/3pe-^^^) + k{q* (3 - P*a) sin(2a;x) + m(a*/3 + {3* a) cos{2ujx) V 
~ Vw(|a|2e2'^^ + |/3|2e-2''^) - K{a*f3 - l3*a) sin(2a;x) - m(a*/3 + (3*a) cos(2cjx) J 

The function mi can be written in the form vi = exp{4i6{x)) where 

e{x) = Arctan ( , ^.'i!. ^ , (/>„ = arga, (/)/3 = arg/3. 

\a; cosh(2K;x + m |a/p|)/ 

It is easy to check that u = v = exp(i/(x)) is an exact solution of (2.1) for any differentiable 
function f{x). It resembles the case of the three- wave equation [48] where one wave of an 
arbitrary shape is an exact solution of the system and the two other waves are identically zero. 
A plot of the solution (3.26) is shown on Fig. 1. The solution (3.26) has a simple spectral 
characterization and an explicitly given analytic fundamental solution of the corresponding 
linear problem. 

Case (ii): /3 = 0, 7/O 

In this case, from (3.13) it follows that the solution can be written in the form 

AiujkQ* ex.p{iujx + ii(K^ — u)'^) + i{(j)a — ^7)} ^ Scjk^ 

{oj — iK)Q'^ ' {uj — iK)Q^' 
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Figure 1. Real and imaginary parts of vi (cf. formula (3.26)) as a function of x. Here k = 1, cj = 10 
a = 1, = 1 + i, 7 0. 




-6 -4 -2 2 4 6 -6 -4 -2 2 4 



Figure 2. Contour plot of Re vi {x, t) (left panel) and Re ui {x, t) (right panel)for a generic soliton solution 
(see formula 3.27) as a function of x and t where a = /3 = 7 = K = tJ = l. 

where 

Q = 2cje'"(^"2"*)+^''l"/^l + (w + iK)e-«(^-2'^*)-i'^l°/^l. 
Case (iii): a = 0,/3/0,7/0 

The solution now can be obtained from the solution in the case (ii), by changing a — )• /3 and 
X — )■ —X. 

In the cases (ii) the solution (3.27) is a soliton of width 1/ n moving with velocity 2uj. The 
corresponding soliton in the case (iii) moves with a velocity —2uj. 

In the generic case, when all three constants are non-zero, the solution (3.13) represents 
a nonlinear interaction of the above described solitons. For k > it is a decay of unstable time 
independent soliton from the case (i) into two solitons, corresponding to the cases (ii) and (iii) 
(see Fig. 2). For k < 0, the solution is a fusion of two colliding solitons into a stationary one. 

Let us now consider the case of imaginary poles, i.e. ^ = in. There exist two essentially 
different cases. 
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-lj>Hr-" : ■ ANf 



Figure 3. Plots of the phase S(x) (left panel) and lmv{x) (right panel) for a stationary doublet soliton 
solution (3.28 ) as a function of x and the phase cf) where k = (Tq = 1- 

1. Suppose po,2 = 7 = 0. Then from (3.23) it follows that |po,i| = IpoM- suffices to pick up 
po,i = 1 and for the third component we have po,3 = exp(z(/9), G M. In this case the 1-soliton 
solution is stationary: 

, , ^ , , /cosh2Kx + sinh2K;rccos(i9 + i(c7o — 2Kxsin w)^ ^ 

u{x) = 0, v[x) = , , r 

\ cosh 2kx + sinh 2kx cos ip — i [ao — 2kx sin 99) 

The function v can be presented as 

f ao — 2Kxsinip \ 

vix) = exp(4«^(x , ^(x) = arctan ; ; . (3.28) 

^ ^ i^v V 7y> \ J Vcosh2Kx + sinh 2ra cos (/^y ^ ^ 

A plot of the solution (3.28) is shown on Fig. 3. In particular, when ip = 0;it, i.e. po,i = ipo,3 
we have for the doublet solution: 



u(x) = 0, v(x) 



g±2KX _|_ ^CTg ^ ^ 



It is evident that the presence of the constant do here is essential since otherwise the solution 
coincides with the vacuum. 

2. Now let us assume po,2 7^ 0. For simplicity we fix po.2 = 1- Then the norms of po,i and po,3 
are interrelated through 

|po,il^ - Ipo.sP = 1- 
This is why it proves to be convenient to parametrize them as follows: 

Po,i = cosh0oe*('^°+'^\ Po,3 = I sinh 6*0 1 e*^'^''-^), ^o,</'o,<^ G M. 
The 1-soliton solution reads: 

m(x, t) = — r^e*('' ^^^°) [sinh(K;x + Oq) cos (p + i sinh(Kx — ^o) sin (p\ , 

A* (A* - 2) 2(2ia - 1) ^iaiia - 1) 

^(^^^)= A2 A + A^ ' 

where 

A(x, t) = cosh^(Kx + ^o) cos^ (p + cosh^(KX — 6*0) sin^ ip — i{cFo — 2K^t + kx sinh 29q sm2(p>), 
a{x, t) = uq — 2K^t + KX sinh 26q sin 2ip. 
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Figure 4. Contour plot oiReu{x,t) (left panel) and Kev{x,t) (right panel) for doublet soliton solution 
of the type (3.29) as a function of x and t. Here k ^ (Tq = 1, Oq = and ~ 0. 



The above solution can be significantly simplified if one assumes that Po,3/po,i £ 
boundary cases the doublet solution reads: 



u{x, t) 
v{x, t) 



sinh(Kx + ^0 



2 (cosh^(Ata; + Op) + ijap - 2kH)) ^^(^2,+^^) 
(cosh2(ra + 6*0) - i{(Tp - 2KH)f 
(cosh^(«;2; + 6*0) + i{cro - 2kH)) [cosh^{K.x + Oq) + iap - 2 - 2iH?t) 



[cosI?{kx + 6*0) - i{cFQ - 2kH)Y 



if Po,3/Po,i > ((^ = 0) and 



sinh(KX — 60 



2 {cosb?{KX - 60) + ijao - 2kH)) ^i(^2i+^^,+^/2) 
(cosh^(Kj; — ^0) ~ ^l^^o ~ 2K?t))'^ 

(cosh^(«;x - 6*0) + i(cro - 2K^t)) (cosh^(KX - 6*0) + iuo - 2 - 2iK^t) 
(cosh^(Kx — ^0) — *(fo — 2K?tyf' 

when po,3/Po,i < ((^ = vr/2). A plot of the solution (3.29) is depicted on Fig. 4. 



u(x, t) 



v{x, t) 



In these 



(3.29) 



3.3 Multisoliton solutions 

As we mentioned above, the dressing procedure can be applied several times consequently. Thus 
after dressing the 1-soliton solution one derives a 2-soliton solution, after dressing the 2-soliton 
solution one obtains a 3-soliton solution and so on. Of course, in doing this one is allowed to 
apply either of dressing factors (3.8) and (3.17). Therefore the multisoliton solution will be 
a certain combination of quadruplet and doublet solitons. 

Another way of derivation the multisoliton solution consists in using a dressing factor with 
multiple poles. For example, if one wants to generate quadruplet solitons one should use 
a dressing factor in the form: 

N 



which is evidently compatible with the reduction condition (3.3). Due to (3.4) its inverse has 
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the form: 



N / 
k=l \ 



Ml{x,t) CMl{x,t)C 



In order to determine the residues of $ one follows the same steps as in the case of a 2-poles 
dressing factor. Firstly, the identity $<I>~^ = U implies that the residues of $ and satisfy 
certain algebraic relations, namely: 



A / mJ cmJc 
^ + X/^M — — + — T"^ 



1=1 



0, k = l,...,N, (3.30) 



to ensure vanishing of the residue of $<I>~^ at A = /i^. Of course due to the Z2 reductions we 
will have an additional set of algebraic relations which are obtained from (3.30) by hermitian 
conjugation. 

As discussed before in order to obtain a nontrivial dressing the residues must be degenerate 
matrices. Thus one introduces the factorization = |nfc)(m/;| and then reduces (3.30) to 
a linear system for \nk) 

I \ V^K I \ H * f {mi\mk) (m,|C|mfc) 

\mk) = }_^Bik\ni), Bik:=fik[- „ , , . C 

^ V W - ^fc W + ^fc 

By solving it one can express the vectors In;), / = 1,...,A through all |mfc) and that way 
determine the whole dressing factor in terms of |mfc). This step contains the main technical 
difficulty in the whole scheme. 

The vectors |mfc) can be found from the natural requirement of vanishing of the poles in the 
differential equations (3.14). The result reads 

\mk{x,t)) = ipo{x,t,fil)\mkfl). 

Thus we have proved that as in the 2-poles case the dressing factor is determined if one knows 
the seed solution ^po{x,t,X). The quadruplet A-soliton solution itself can be derived through 
the following formula 



a + Y,iMk + CMfcC) Lf ) U + J^(Mfe + CMfcC) . 

V k=l / \ k=l / 



From all said above it follows that the algorithm for obtaining the A^-soliton solution can be 
presented symbolically as follows 

^ (|mi), . . . , |m,v)) ^ (|ni), . . . , |n^)) ^ (Mi, . . . , M^) ^ L^^\ 

Similarly, one is able to generate a sequence of N doublet solitons by making use of the 
following factor 

N 



^^ . \^.fPlix,t) CPl{x,t)C\ 



Following the same steps as in the 1-soliton case one can convince himself that the vectors |q; 
involved in the decomposition Pi = |qz)(pi| satisfy the linear system: 



N 



s=l 
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where 

T>ss ■= t(Ts C., T^sl ■= m ; , s 7^ t- 

2 \ Ks - Kl Ks + Kl J 

On the other hand the vectors \pi) and the real functions u/ are expressed in terms of the initial 
solution ■00 and its derivatives with respect to A as follows: 

\pi{x,t)) = ipo{x,t,-iKi)\pifl), 

ai{x, t) = -Kiipifili/jQ^ix, t, iKi)ipo{x, t, iKi)C\pifi) + ai^o, 

where \pifi) and ai^ for / = 1, . . . , are constants of integration to parametrize the solution. 
The components of |p;) are not independent but satisfy 

(p,|C|p,)=0 (pz,o|C|pi,o) =0. 

In order to derive doublet A-soliton solution one uses: 

^S'^ = (^a + J2{Pi + cPiC?j Lf^ + Y^iPi + CPiC) 

In this case the diagram which describes the algorithm to obtain the multisoliton solution reads 
lS°) ^ {\pi),ai}l, ^ {|q,)}ili ^ {Pi}l, ^ L«. 

It is clear that apart of the 'pure' multisolitons discussed above there is a 'mixed' type of 
multisolitons, i.e. multisolitons to contain both species of 1-soliton solutions. In order to obtain 
such solution one needs to use a dressing factor possessing both types of poles. 




4 Generalized Fourier transform 

In this section we aim to develop the generalized Fourier transform interpretation of the inverse 
scattering method for the Lax operator (2.3). For this to be done we are going to use some ba- 
sic notions like Wronskian relations, 'squared solutions', recursion operators etc. The 'squared 
solutions' also known as adjoint solutions are generalizations of exponential functions in usual 
Fourier analysis, namely they are a complete system. They occur naturally in Wronskian rela- 
tions which interrelate the scattering data and the potential Li(x) (as well as their variations) 
in a form which resembles a scalar product and is called skew-scalar product. The Wronskian 
relations allows one to express the expansion coefficients of Li{x) and its variation along these 
generalized exponents in terms of the scattering data and its variations as we shall see later 
on. Within this framework the role of the operator id/dx whose eigenfunctions are exp(ifcx) is 
played by the recursion operators. We shall present two substantially different ways of obtaining 
the recursion operators and compare them. 



4.1 Recursion operators and integrable hierarchy 

Recursion operators appear in many situations in the theory of integrable systems [1, 20, 41]. 
Their existence is tightly related to the fact that each integrable NLEE belongs to a whole family 
of infinite number of integrable NLEEs associated to the same Lax operator L. One can also 
assign to this family a whole infinite set of symmetries and integrals of motion (Hamiltonians). 
In this subsection we are to derive the recursion operators by analyzing the interrelations between 
members of the hierarchy of NLEEs. In doing this we are following the ideas in [24, 46]. 
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The hierarchy of integrable NLEEs associated to L is generated through the zero curvature 
representation which we rewrite in the original Lax form 

iLt = [L,V]. 

Each choice of V corresponds to an individual member of the hierarchy. Following the ideas 
in [24] we interrelate two adjacent flows V and V through the equality 

V{X) = k{\)V{\) + B{\), (4.1) 

where k[\) is a scalar function and the remainder B has the structure of the Lax operator A{\) 

B{\) = XBi + X^B2. 

All quantities above must be invariant with respect to the action of reductions (2.7), (2.8). This 
is why we pick up k{\) = while the operator B obey the condition 

B\X*) = B{X), CB{-X)C = B{X). 

This means that the coefficients of B are hermitian matrices to fulfill Bi G s[^(3), B2 G sf{3). 
After substituting (4.1) into the Lax representation one obtains the following relation 

iLr = iX'^Lt + [L,B], (4.2) 

where r is the evolution parameter of the flow V. Then the recursion operator TZ is defined as 
the mapping 

Li^T = TiLi^f 

In order to find it we need to solve (4.2) and thus calculate the remainder B. 

Since relation (4.2) holds identically with respect to A it splits into the following set of 
recurrence relations: 

iLi,t + [Li,52] =0, (4.3) 

zB2,, + [Li,5i] = 0, (4.4) 

Li,r - Si,^ = 0. (4.5) 

In accordance with our conventions explained in the Preliminaries (see formulae (2.10)-(2.12)) 
the coefficients Bi and B2 have splittings 

Bi = Bi + hLi, {Bi,Li)=0, 

B2 = Bi + 62^2, {Bi.L2) = 0, (4.6) 

of terms B^ and B2 not commuting with Li and their orthogonal complements. After substi- 
tuting (4.6) into (4.3) it becomes possible to express B2 through Li^t as follows: 

Bi = -i8.d-^lLi^t, (4.7) 

where ad^^ acts on the orthogonal parts only, i.e. the kernel of ad/,^ is factorized. In this case 
it can be proven the presentation below holds true 
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The coefficient 62 is determined from relation (4.4). Indeed, after plugging (4.6) into (4.4) 
we get 



idccB^ + ib2,xL2 + ib2L2,x + [Li,Bi] = 0. 



The Li-commuting part of (4.8) is separated by taking ( • , L2) in both hand sides of the equation 
to give 



-,b2,x + {dxBt,L2) =0 



-d-\d,B^,L2), 



(4.9) 



where we have denoted by ^ the integral dy. On the other hand by projecting (4.4) with 
vr = ad^^^ ad-Li one extracts its orthogonal part 

ind^B^ + ib2L2,x + [Li,Bi] = 0. 
Taking into account (4.9) one can find out the following recurrence relation 



Bi = A2Bi, A2 = -iad^i (^TTd^i-) - h2,xd-\d^{-),L2)^ . 



(4.10) 



Formally speaking, one should write instead of A2 for the integro-differential operator in- 
troduced above since the lower integration limit in differs. This choice of signs, however, 
is unessential for the considerations in the current subsection and this is why prefer a more 
simplified notation. 

After a similar treatment of (4.5) we obtain 



26i,a. + (5,5f,Li) =0 



bi = -\d-'{d,Bi,L,), 



L,^r = i^dL,AiBi, Ai = -iad^i (^7rd,{-) - Li^.^d-' {d^{-), Li)^ . (4.11) 



Finally combining (4.7), (4.10) and (4.11) one derives 

Li^T- = ad^i A1A2 adj^^^ Li^t ^ 7^ = ad^^ A1A2 ad^^^ . 



(4.12) 



The recursion operator can be also presented conveniently as a 4 x 4 matrix to map the column 
vector {u,v,u* ,v*)J to the column vector {u,v,u*,v* )l. Then it can be verified that TZ obeys 
the splitting: 

7^ = 7r^(i3loc + Snonl) + /C. 

The factors vr and A look as follows 



vr 



/ 1- 


-|n|V2 




-uv*/2 


-u^l2 


-uv/2 \ 




-u*v/2 


1 


-\v\^/2 


-uv/2 


-vy2 




{u*Y/2 




-u*v*/2 1 - |n|V2 


-u*v/2 


v - 


u*v*/2 




-(i;*)V2 


-uv*/2 


1 - |f|V2 / 




( ° 


— V 


-n \ 






idx 


—u 





V 


dx- 




* 

V 





u* 






\ 


u* 


-V* ) 
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The matrices 



loc 



/ -4u* 2i;(3|t;|2 - 1) -6u*v'^ \ 

2t;*(l - 3|up) 6u{v*y 4u 

-4u* 3u(|w|2 - |Mp) -v{l + 6\u\'^) 

\ -u* -V* 



J 



and 



B 



3z 



nonl 



/ 2u*v \ 
2uv* 
lup — \v\^ 

\ -1/3 / 



are a local and a nonlocal part respectively originating from the orthogonal part B2 and the 
Li-commuting part ^2-^2- Finally 



1 



d-'[{u\v\u,v)A{Bxoc + B, 



nonl J 



V 
U* 

\v*l 

is connected with the term h\L\ in the splitting of B\. 



4.2 Wronskian relations and 'squared solutions' 

Wronskian relations [5, 6] provide an important tool for analyzing the relevant class of NLEE 
and the mapping 5^ : 2)t — >• T, where OJt is the set of allowed potentials of L (in our case L\) 
and T is the minimal set of scattering data. 

In deriving them we will need along with the first equation in (2.16) also two other related 
equations: 



' dx 
.d&x 

' dx 



- Ax(x, A)Li(x) = 0, x{x,\) = x ^{x,X), 
+ XLi{x)6x{x, A) + \5Li{x)x{x, A) = 0, 



where the variation of xi^-, A) is due to the variation 5Li{x). 

The first type of Wronskian relations interrelates the asymptotics of FAS with Liand its 
powers as shown in the examples below 



(4.13) 
(4.14) 
(4.15) 



J —00 

/oo 
dx xLi,xX{x,X), 
-00 

/•oo 

X^2(x)x(a;,A)r = / dxxL2,xX{x,X)- 



We remind the reader that J = diag(l,0, — 1) and I = diag(l/3, — 2/3, 1/3) are the diagonal 
forms of (the asymptotics of) Li and L2 = -^1 — 21/3 respectively, that we introduced at the 
very beginning in section Preliminaries (see formulae (2.13), (2.14)). 

A second class of Wronskian relations connects the variation 5xix, A) with the variation of Li, 

i.e. 



X^x{x,X)\ 



00 
—00 



iX / dx x6Lix{x, X). 
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In the left hand sides of the Wronskian relations are involved the scattering data and its variation 
while the right hand sides can be viewed as Fourier type integrals. To make this obvious let us 
take the Killing form of the Wronskian relation (4.13) with a Cartan-Weyl generator Ea and 
use the invariance of the Killing form 



oo 



i{xJxix,X) - J,Ea) = A / dx ([Li, J], 60(3;, A)) . 

J —00 

The quantities 

ea{x,X) = xEaX{x,X) 

introduced above are called 'squared solutions'. Due to the fact that we have two FAS x^ix, A) 
and x^{x,X) we obtain two types of 'squared solutions' e^(a;,A). Similarly, taking the Killing 
form in (4.14) and (4.15) we find 



{xLix{x,X),Ea)\_^= dx {Li^:r,ea{x,X)) , (4.16) 

J ~oo 

/oo 

dx {L2,x,ea{x,X)) , (4.17) 
'OO 

;>oo 

ixSxix, A), Ea) r^ = iX dx {6Luea{x, A)) . (4.18) 



00 



We are interested more specifically in variations that are due to the time evolution of Li(x), i.e. 

dL 

SLi{x) = Li{x,t + 5t) - Li{x,t) ~ ^t-Q^- 
Therefore up to first order terms of 5t we obtain 



{xXt{x, X),Ea) = iX I dx {Li^t, ea{x, A)) . 



Now we can explain why the Wronskian relations are important for analyzing the mapping 
5" : 9JTl^ — )• T. Indeed, taking x{x, A) to be a fundamental analytic solution of L we can express 
the left hand sides of (4.16) and (4.17) (resp. (4.18)) through the Gauss factors S^, and 
(resp. through the Gauss factors and their variations). The right hand side of (4.16) and (4.17) 
(resp. (4.18)) can be interpreted as a Fourier-like transformation of the potential Li{x) (resp. 
of the variation 6Li{x)). As a natural generalization of the usual exponents there appear the 
'squared solutions'. The 'squared solutions' are analytic functions of A. This fact underlies the 
proof of their completeness in the space of allowed potentials, as we shall see. 



4.3 The skew-scalar product and the mapping VJIl^^ — >• %j 

It is obvious, that only some of the matrix elements of the squared solutions ea{x, A) contribute 
to the right-hand sides of the Wronskian relations. To make this clearer we will use the Z2- 
grading of the Lie algebra which hints that we should split the squared solutions as in (2.10), 
namely: 

ea{x,X) =naix,X)+ICaix,X), ^„ (x. A) G 5[° (3) , /C^ (x. A) G sl^ (3) . (4.19) 

In addition, according to (2.11), (2.12) each components T-La{x, A) and /Cq,(x, A) can be split into 

3 

Ua{x,X) = Ha{x,X) + haL2{x), ha = - (-L2 (x) , ^q, (x. A)) , 
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}Ca{x,X) = Ka{x,X) +kaLi{x), ka = - {Li{x) , ICa{x, X)) , (4.20) 

where Ha{x, A) and Ka{x, A) do not commute with Li(x). It is not difficult to realize that only 
Ha{x, A) and Ka{x, A) contribute to the right-hand sides of the Wronskian relations. 
In what follows we will make use also of the skew-scalar product 

/oo 
dy{Xiy),[L,iy),Yiy)]), 
-oo 

where X and Y are functions with values in s[(3)/keradL^ and vanishing for x — )• ±00. 

We will denote the linear space of such functions by TIli- Note, that we can express the 
right hand sides of the Wronskian relations using the skew-scalar product: 

ixLixix, X),Ea) = [[ea{x, A), ad^^^ Li^^:]] , 

{xL2X{x, X),Ea) |!^_oo = ha{x, A), ad^^^ L2,x]] , 
{xSx{x,\),Ea) |!^„oo = iX\[ea{x, X), adll 6Li]\, 
{xXt{x,X),Ea) |^_oo = iX[[ea{x,X),adll Li^t]]- 

We should also point out that the skew-scalar product is non-degenerate on TIli- 

We finish this subsection by listing the effective results from the Wronskian relations. Below 
a e A+: 

pW'±(A) = {D^f^JT^D^,E±a) = -[[^d-^lLi,x,Ki^{x,X)]], 

pi2)-±(A) EE {D^f^IT^D^,E^^) = -[[^dllL2,.,Hi^{x,X)]]; (4.21) 
rW'±(A) = {S^JS^,E^^) = [[ad^iLi,,,/^±„(x,A)]], 

r(2).±(A) ^ {S^IS^,E^^) = \ladllL2,.,H$^{x,X)]]; (4.22) 
S'ptiX) = {D^f^6T^D^,E±^) = -iX[[s.dll6Li,Ki^{x,X)]\, 

5'T^iX) = {S^6S^,E^o.)=iX\l^dll6LuK$^ix,X)]]. (4.23) 

In deriving the above formulae, besides the splitting (4.19) we also used the fact that 
ad^^^ 6Li,adll Li,^ G 5[°(3) and ad^^^ L2,x G sl^(3). 

The quantities p^^'^ , Ta^'^ , k = 1,2, can be viewed as reflection coefficients. These formulae 
provide the basis in the analysis of the mapping Tlii Tfc, A; = 1, 2. Indeed Li,^ G ^Li, while 
the sets of coefficients 

Ti = 1? U \ lP EE {t(^)'±(A), a G M, a > 0}, 
T2 = U = {pi^^'^W, A G M, a > 0} 

completed with additional sets characterizing the discrete spectrum of L are candidates for the 
minimal sets of scattering data of the Lax operator. 

The above considerations show, that following [1] we can treat the 'squared solutions' as 
generalized exponentials. Therefore the mapping 971^, can be viewed as a generalized Fourier 
transform. 

Of course, the justifications of the above statements must be based on a proof of the com- 
pleteness relation for the 'squared solutions' which will be presented in the next section. 
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4.4 Recursion operators — an alternative approach 

Our goal in this subsection is to present an alternative definition of the recursion operators. 
We introduce them as operators whose eigenfunctions are the 'squared solutions'. More 
precisely, we shall see that the following equality holds: 

A±i^^(x,A) = A2K±(x,A). 

Our consideration starts from analysis of the equation 

idxCa + X[Li{x),ea{x, A)] = 

satisfied by each of the 'squared solutions'. We have omitted the superscript it in the notation of 
the squared solutions since it does not play a significant role in most of our further considerations. 
We shall use it again at the end of subsection when it will become essential. 

Due to the grading condition (2.9) the components Tia and /C^ are interrelated through the 
system 

z9x^a + A[Li(x),/C„(x,A)] =0, idxlCa + X[Li{x),na{x,X)]=0. (4.24) 

After we extract terms proportional to Li from the above equations we obtain 

3 

{L2,dxHa) + d^K = ^ K = hafl - -d~^{L2,dxHa), 
{Li.dxKa) + dxka = =^> ka = kafi - ^d^^ {Li, d^Ka) , 

where hafl and ka,o are some integration constants. On the other hand the orthogonal part 
of (4.24) reads: 

indxHa + ihaL2,x = -X[Li{x), Ka], iirdxKa + ikaLi^x = -X[Li{x), Ha\. 

After substituting and k^ in the equations above we get 

kiK^ = XHa - ka,oA2L2, (4.25) 
K2H^ = XKo,-KflKiLi, (4.26) 

where the operators A2 and Ai are given in (4.10) and (4.11) respectively. This is the right 
place to mention that remarkably the operators Ai^2 are invertible. Indeed, they admit the 
factorization 

kiX = -is.dllJidxX, JiX = X- hi,xd-\L^,X), 

A2Y = -iadll ^2a.y, J2Y = X- ^L2,xd~\L2,X). 
It is not hard to see that the operators and ^2 can be inverted as follows 

Jf^X = X + hi^xdxHLuX), J^-'Y = X + h2,xdxHL2,X). 
Therefore Ai and A2 are invertible as well, namely 

A^^A: = idx^J^^ adL, X, A^^X = id'^J.^^ ad^, Y. 

Let us apply A2 to (4.25) and Ai to (4.26). After restoring the index notation where needed 
the result reads: 

A^AfKt = X^Kt - A/ia,oAf Li - fc„,oAf Af L2, 
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This is the right place to remind the reader that the constants hafl and kafl are determined 
by the asymptotic of the relevant 'squared solution' for x — )■ oo (or for x — >• — oo), depending on 
the proper definition of the recursion operator. More detailed analysis shows that for each A^ 
and A^ there exist certain roots for which the constants /ca,o a-nd hafl vanish. More specifically 
the following equalities holds: 

A+i^|jx, A) = Ai/|„(x,A), Ari^|jx,A) = AF|Jx,A), 
A+F|,(x, A) = AK|„(x, A), A2-i^±„(x, A) = AK|,(x, A) 

for all roots a > and therefore: 

A+A+K±„(x, A) = X^K^^ix, A), A2 Ari^|,(x, A) = X^K^^{x, A), 
A+A+F±,(x, A) = X^H^^ix, A), ArA2 //|„(x, A) = X^H^^{x, A). 

The operator A^ introduced as 

K^X = KfK^X, Xgs[°(3), Ay = AfAfy, yGsl^(3) (4.27) 

is the recursion operator we have been looking for. If we compare (4.27) and (4.12) we can 
conclude that both approaches lead to compatible results which coincide up to some conjugation 
by adii • 

5 Spectral theory of the recursion operators 
5.1 Asymptotics of the fundamental analytic solution 

For performing the contour integration, in order to derive the completeness relations for the 
squared solutions, one needs the asymptotics of the squared solutions for | A| — )• 00. The starting 
point here is the asymptotic expansion of the FAS for A — >• 00: 

00 

X^s{x,X) = t + ^Xk{x)X-^ . 

k=l 

Introduce: 

Xas(a;, A) = g-\x)xas{x, A)e^-^^^, 
which satisfy: 
■ dx 



dx 



+ ig ^gxXasix, A) + A[J, Xas{x, A)] = 0. (5.1) 



Together with properly chosen asymptotic conditions Xa.s{x,X) the last equation will provide 
the asymptotics of the fundamental analytic solution. Therefore it will allow an asymptotic 
expansion of the form: 



Xas(2;, A) = ^ A "xkM^)- 

s=0 



This will lead to a recurrent relations for the expansion coefficients Xk,asix)- Inserting this 
expansion into equation (5.1) for the first two coefficients of Xasix, A) we get: 

[J,XO,as(x)] =0, (5.2) 
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^-^p - iQx ^9Xo,as(a;) + [J, xi,as(a;)] = 0, (5.3) 

.C^Xl.as . _i - 



From (5.2) it follows that Xo,as{x) must be a diagonal matrix and using the diagonal part of (5.3) 
one can derive it in an explicit form: 

Xo,as(x) = diag(e^(-),e-2^(-),e^(-)), p{x) = \ T dy{u*Uy + v*Vy), (5.4) 



iboo 



and for the off-diagonal part of Xi ^3(2^) we have (see equation (2.15)): 

V2{uvx - v'Ux)e-^(^') {u*u^ + 7;*t;^)e''(^) /2 

x{ as(^) = ^ ( V2iu*v* - z;*OeP(^) -V2{u*v* - v*ul)e'^p'^^'> 



2 



As a result the asymptotic behavior of xi^^ ^) for A — )• 00 is given by: 

x{x,X) ~ <7(2;) (xo,as(a;) + YXi,as(x) + --- je-^-'^^ (5.5) 

A— !>oo y A / 

Note, that the fundamental analytic solution x^ does not allow a canonical normalization for 
|A| — )• 00. This difficulty can be overcome by applying a suitable gauge transformation. 

Using the asymptotic behavior (5.5) of the FAS, one can derive the asymptotics of the squared 
solutions ea{x,X) = {x{x, X)EaX~^{x, X)Y . Skipping the details: 



etix, A) e^('''"^^-g{x)x.s{x)E^x^Hx)9-\x) 

A— >oo 

= e*W")^^-+3^W(/,")^(2;) (^Ea + 0{X-^)) g-Hx)- 
5.2 Completeness of squared solutions 

Below we will derive the completeness relation for the 'squared solutions' for a class of poten- 
tials Li{x) of the Lax operator L which for every fixed value of t satisfy the following conditions: 

CI) Li{x)—Li^± is complex valued function of Schwartz type, i.e. it is infinitely smooth function 
of x falling off for |x| — )• 00 faster than any power of x; 

C2) Li{x) is such that the corresponding functions mf{X) and m^(A) have at most finite num- 
ber of simple zeroes. For simplicity we assume also that mf{X) and m^(A) (resp. (X) 
and 1712 (X)) have no common zeroes; 

C3) the corresponding functions m^(A) and m^{X) have no zeroes on the real axis of the 
complex A-plane. 

Remark 2. It is well known that the zeroes of the functions mf{X) and m^(A) determine the 
discrete spectrum of L. The proof of this fact comes out of the scope of the present paper. For 
other classes of Lax operators it is well known, see e.g. [11, 20]. 

Below for definiteness we introduce notations for the discrete eigenvalues. As we already 
mentioned, due to the two Z2-reductions, they are of two types. Let there be A^i discrete 
eigenvalues of generic type which come in quadruplets (see Fig. 5). We will denote by A^, A^, . . . , 
A^^ those of them that are in the first quadrant, i.e. < argA^ < vr/2, 1 < A; < A'^i. The 
eigenvalues A]^^_,_^ = — (A^)* belong also to C+ but are in the second quadrant. Besides we have 
also XJ = (A^)*, j = 1, . . . , 2Ni which lie in C_. The second type of discrete eigenvalues are 
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Figure 5. The bold line defines the continuous spectrum of L, the regular lines show the integration 
contours 7±; by 7±_oo we denote the 'infinite' semi-circles. By o we denote the discrete eigenvalues 
relevant for one generic soliton which come in quadruplets ±A^; the second type of solitons correspond 
to purely imaginary pairs of discrete eigenvalues denoted by x . 



purely imaginary, i.e. \2Ni+a ~ "^^a, a = 1, . . . , A''2; obviously A^^^_,_^ = —iva- Therefore the Lax 
operator L will have 4A^i + 2N2 discrete eigenvalues. 
Consider the Green function: 

G^{x,y,\) = Gt{x,y,\)e{y - x) - Gt{x,y,\)e{x - y), (5.6) 
Gf{x,y,\)= ^ e%^{x,\)®e%^{y,\), (5.7) 

aGA+ 

r 

Gf(x,y,A)= e±„(x,A)®e^^(y,A) + ^/i^±(x,A)0/i^±(y,A), (5.8) 

aSA- j=l 

hf{x,\) = x^{xA)HjX^{x,\). (5.9) 

In order to derive the completeness relations for the squared solutions, one needs to apply 
a contour integration to the integral: 

■^G{x,y) = ^^j> dXG^{x,y,X) - j) dXG~{x,y,X), 

where the integration contours 7-1- = R U 7±,oo are shown on Fig. 5. 

Let us recall that the Lax operator L has two types of discrete eigenvalues whose number 
is 4A'^i + 2N2; the first type correspond to generic solitons and come in quadruplets iA^; the 
second type correspond to purely imaginary pairs of discrete eigenvalues ±ii^a- Thus according 
to Cauchy's residuum theorem 

JGix,y) = V Res G+{x,y,X) + Res G-{x,y,X) , 

presuming that the radii of the 'infinite' semicircles are large enough so that all discrete eigen- 
values are inside the contours 7-1- . Condition C3 above ensures that G^{x, y, A) has no poles on 
the continuous spectrum. 
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Next we integrate along the contours: 
1 

^G{x,y) = — J dA(G+(x,y,A)-G-(x,2/,A)) 

In order to calculate the integrals over ^±,c)o it is enough to use the asymptotic behavior of 
the Green function for A — )• oo. From equations (5.4) and (5.5) we have: 



Q± (^x,t,X) ~ ^i{J,a)X{x-y)+3(pix)-p{y)){I,a) 

X {g{x) ® g{y)) {E±^ ® E^^ + 0{\-^)) {g-\x) g-\y)) 

and 

r 

G%^{x, t,X) ^ Y.^9{x) ® g{y)) {Hj HJ + ©(A'^)) {g-\x) g-\y)) 

^ ^i{J,a)X{x-y)+3{p{x)-p{y)){I,a)^g^^^ ^ ^^y^^ ^^^^ ^ ^ 0{\-^)) {g{x) (g) g{y)). 

oeAouA^ 
As a result we get: 

^ / dA Gi(x, y^^)-^-f dX G-,{x, y, A) 

"''7+, CO "''7— ,oc 

= 5(x-y)(5(x)0 5(y)) ^ (^«(»^_^-^_a®^a)(5(a;)®5(y))- 

The next step is to simplify the integral along the continuous spectrum; namely we will 
calculate the jump of the Green function with the result: 

G^{x,y,X) - G~{x,y,X) = G^{x,y,X) - G^{x,y,X) 

= Yl (e^(x,A)®et^(x,A)-eI„(x,A)®e-(x,A)). 

The proof of this fact goes as follows. From the expression for the Green function (equa- 
tions (5.6)-(5.9)) and using the relation 9{x — y) + 0{y — x) = 1 we easily get: 

G+(x, 2/, A) - G-(x, y, A) = G+(x, y, A) - Gr(x, y, A) 

- e{x - y)(G+(x, y, A) + G+(x, y, A) - ^^(x, y, A) - G2 (x, y, A)), 

so it will be enough to prove that the coefficient in front of 9{x — y) vanishes. Indeed, we have: 

Gt{x,y,X) + G^{x,y,X) - G^{x,y,X) - G2 (x,y, A) 
= (x+(x, A) x+(y, A))n2(x+(x, A) ® x+(y, A)) 
- (x"(x, A) ® x'(y, X))ii2{x'{x, X) x~{y, A)). 

Here 

2 

n2 = V (^^a $5 ^-a + C$ ^,) + V 
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and Hi = J, H2 = J4. Note that 112 is the second Casimir endomorphism for the algebra s/(3). 
It remains to use equation (2.25) and the basic property of 112 that states: 

(G(A) G(A))n2(G'(A) G{X)) = U2, 

for any element G{X) of the corresponding group SL{3). 

We evaluate the contribution from the discrete spectrum assuming that the squared solu- 
tions e^(x, A) have at most simple poles at the points of the discrete spectrum. Therefore in 
the vicinity of the eigenvalues A^ we use the following expansions of the 'squared solutions': 

e^(x,A) = (A-A^r + 
where pa are integers and pa = —p-a- Thus we obtain: 

Res^G+(x,y,A)= ^ ^ e^^,^ + ® e±^^,). 

Therefore the completeness relations takes the form: 
S{x - y) {g{x) ® g{y)) Ho {g{x) ® g{y)) 

/oo 



2^ 



2N1+N2 

A:=l aeA+ 

+ e:„ ® e- ,{y) + el^ e- ,{y)) , (5.10) 



where 



no = 2^ — 



(a, J) 



and 



5.3 Expansion over the 'squared solutions' 

Using the completeness relations of the squared solutions, one can expand any generic ele- 
ment F{x) of the phase space Ai over the complete sets of squared solutions. We remind 
that F{x) is a generic element of A4 if it tends fast enough to zero for |x| — )• 00 and is 'perpen- 
dicular', i.e. 

(F(x),Li) = (F(x),L2) =0. 

It can be parametrized as follows: 

F{x)=g{x) Yl Fa{x)E^g{x). 
oeA+uA_ 
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Let us now multiply both sides of the completeness relation (5.10) by 1 [Li{y),F{y)] on the 
right, take the Killing form of the second factors in the tensor product and integrate over y. In 
the left hand side after some manipulations we get: 

dy {(gix) g{y))Uo{g{x) g{y)), U ^ [Li{y) , F {y)]) ^ 5{x - y) 
= ^(^) E r\\3{x) dy{g{y)E.M,[Li{y),F{y)])6{x-y) 
= gix) Yl Fa{x)E^g{x) = F{x). 



ieA+uA_ 



The same operations applied to the right hand sides of (5.10) replaces each of the terms e^(x, A)(8' 
e^Q,(y, A) by e^(x, A) [[e^(?/, A), -F(y)]] • Thus the expansion coefficients of F{x) are provided by 
its skew-scalar products with the 'squared solutions'. 
The result is: 



nx) = ^ y^'iA ^ {etix,Xh+J\)-eZ^{x,Xhp.JX)) 



oeA^ 
2N1+N2 

+ E E (5-11) 



where 



7F;a(A) = [[e^,(y,A),F(y)]], ^|,«;fc(x) = e±„^,(x)7|,,^, + e±^^,(x)7^,„^„ 

l$,a;k = Hady)^^(y)l' '^F,a;k = ¥Ta;kiy) ^ Hv)]] ■ (5-12) 

Similarly, exchanging x and y, multiplying both sides of the completeness relation (5.10) by 
[Li{y), F{y)] U on the right, take the Killing form of the first factors in the tensor product 
and integrate over y we get: 



P(^) = -^ J J^T. (4(x,A)^+„(A)-e:„(x,A)^^^„(A) 



asA^ 
2N1+N2 

E E(^F^.,.(^) + ^F-«,fc(^))' (5-13) 



where 



= [[eLiy, A), F{y)]] , y±„^,(x) = e^^.^,{x)9%^.^, + 4,.^,{x)e%^.^„ 

&%a;k = hL;kiy)^F(y)]], e%^.^, = [[et^,,{y),F{y)]]. (5.14) 
The completeness relation (5.11) allows one to prove the following 
Corollary 1. 

i) The function F{x) = if and only if all its expansion coefficients (5.12) vanish, i.e.: 

lF;-aW = lF;aW = « ^ A+; 

ZF,a;ki^) = = k = l,...,N, (5.15) 

where N = 2Ni+ N2. 
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ii) The function F{x) = if and only if all its expansion coefficients (5.14) vanish, i.e.: 

(^t;-aW = <^F;aW = 0, « G A+; 



A; = 1, 



iV. 



Proof. Inserting F{x) = into the r.h.s. of the inversion formulae (5.12) we obtain the rela- 
tions (5.15). Next, let us assume that equation (5.15) holds and let us insert it into the r.h.s. of 
equation (5.11). This immediately gives F{x) = 0. Thus i) is proved, ii) is proved analogously 
using equations (5.14) and (5.13). ■ 

In other words we established an one-to-one correspondence between the element F{x) G M 
and its expansion coefficients. 

Now, let us take F{x) = adJ^^Lj^xix), j = 1,2. Then the corresponding expansion coefficients 
are given by: 

£^'^{X) = [[ei^{x,X),s.dll Lj4, 
ri^'^'^(A) = -[[e^«(x,A),ad^iL,,.]], 

From the Wronskian relations we have: 
pi2)-±(A) = 3(i)±fT/rTZ)±,i?^„), 



o^^^'^ - 



'a;k 



Taking into account that ad^^^^ Li^^i^^) G sl^^^ and ad^^ L2^x{x) G st*-^^ we obtain the following 
expansions of these functions over the 'squared solutions' (for details see the appendix): 



ad^^^ Li^x{x) 



2^ 



ad^^^ Li^x{x) 



/oo 
dX {a,J){s+{X)H+ix,X) + s-iX)HZ^ix,X)) 

N 

/oo 
dX ia,J){r+iX)H+^ix,X)+r-{X)H-{x,X)) 

N 

+ E E («'^)«.^-a;^-(x)-r;;,^-;fc(x)). 



k=l oeA 

2^ 



k=l aeAa 



Similarly: 



I 

2^ 



/oo 
E (r('^+(A)i^^(x,A)-r(2).-(A)Kl„(x,A)) 



N 



-EE (rS^Kk(-) + r!S-KZ^^k(-))^ 

k=l aeA+ 
■ /-co 

ad^iL2,.(x) = -- J JX Y (pP+(A)i^^Jx,A)-p(2),-(A)i^^-(^,A)) 



2tt 

N 



leA 



+E E (pSr^^;.(-)+<ri^ 

k=l aeA+ 



a:k 
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Here 



±(a I)s^(X) 


for rv 




^f^=Res4^)'-(A) 


±(a,I)r±(A) 


for a G (5^, 









and Sf = {ai, Q2}. 

Next, taking F{x) = ad^^ SLi{x) we derive the corresponding expansions of the variation 
of Li{x) over the squared solutions. Thus we obtain the following expansions for ad^^^ SLi{x): 

1 f'°° rl\ 

i^ll5L^{x) = -- -Y. (//+(x,A)5'r(i)+(A) -ifrJx,A)(5r(i)-(A)) 



27r 

N 



EE U 

k=l oeA+ I k 



+ 



1 r°° fix 

l^ll5L,{x) = -j ^- J2 {hUx,\)S'p^^^+{X)-H~{x,\)Sp2^-{\)) 



+ H_^.k{x)T^.f^5X^ 



aeA+ 



N 



+ EE U 

k=l aeA+ I k 
1 



+ 



X: 



PaJK 



+ K:ki^)Pa-JK 



The corresponding expansion coefficients arc given by: 



6't^{X) = SstiX) for a G <5+ <5'r± (A) = 6st, ± ^(^44 - 



S'ptW = SrtiX) for a e <5+ ^^^(A) = Jr^ ± -(<5r± r± - ). 

The calculation of the expansion coefficient for the discrete spectrum is explained in the ap- 
pendix. 

Presuming that the variation of the potential is due to the time evolution of Li(x): 

SLi{x) = Li{x,t + St) - Li{x,t) ~ St^, 

one can get similar expansions for the time-derivative ad^^^ Li^t{x): 

1 f'°° fix 

iadZ,iLi,i(x) = -y_^- X; {Str+{X)H_^{x,X)-StT-{X)H^{x,X)) 



EE^^ 

k=l aeA+ 



Polynomial Bundles and GFT for Integrable Equations on A. Ill-type Symmetric Spaces 33 



+ ^-a;fc(^)^a;fcV 



and 



1 1'°° d\ 

ad^iLi,i(x) = - j ^^Y. {StptiX)H+^{x,X) - Stp-{X)H-ix,X)) 



TV 



k=l aeA+ 
1 



Pa;kK,t \ 
XI ] 



xt 



+ ^^.^{x)p^.^X 



k,t 



where the expansion coefficients 6tTa^'^ and StPa'^'^ are given by: 



StptW 



foraG5+, 5tT±(A) = 

drt{X) 



dt 



dt 



for a G (5+, ^tptzW 



03 



dt 





\ 


rdst. 








V dt 


''02 


± 


\ 




r± 






[ dt 


02 



dt 
"1 dt 



The details of calculating the expansion coefficients relevant for the discrete spectrum of L are 
given in Appendix A. 

The completeness relation (5.10) and the expansions that we constructed above can be viewed 
as spectral decompositions for the recursion operators and (4.27). Indeed, one may check 
that 



and 



A+Hi„{x,X) = X^Hi^ix,X), 



K+K^^{x,X) = X^K^^{x,X), 

^^K^a;ki^) = {X^,fKi^,,{x), 



A-H^,ix,X) = X'H^^{x,X), 



^'-HL;ki^) 
^-H^a;k(^) 



iKrH^a;ki^)+2X^H^^.^, 

iKfH^.;ki-) + 2XtH^^.^,{x) 



(5.16) 



A-i^±,(x,A) = A2K±Jx,A), 



~^^KlAx) 



i&K^a:ki^) + 2XfK 



In other words we have proved the completeness of the set of eigenfunctions of A± and A± (4.27). 



5.4 Conjugation properties of the recursion operator 

Here we derive the recursion operators, that are adjoint with respect to the skew-scalar product, 
i.e. we define A^ 2 by the relations: 

[[X2,AfX^]] = [[{AfYX2,X^]], [[Y2,AfY,]] = [[ (A±)* F^, I'l]] , (5.17) 

where Ai,2 = X^^ e s[^(3) and ^1,2 = G s[°(3). 

In doing this we use several times integration by parts and the following properties of the 
operator ad^^^ 

[Li,ad^i Ai,2] = Ai,2, (Ai,ad^i A2) = -(ad^^^ Ai, A2), 
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[Li, ad^i yi,2] = 11,2, ad^i Y2) = -{adll Yi, 

The derivation of (A^)* goes as follows: 

[[X2,AfX^]]=- 



2 ■ 



dx ( X2 



Li,zad^^ — 



+ 



dx (X2 



dx ( ^ ^ , Xi 



—00 
00 



da; ( ad 



dx '^V"^ 
1 8X2 



' ±00 



dy ( Li, 



±00 



dXi 
dy 

dy {Xi,Li^y) 



,[^1,^1] 



3X2 



Li 

-l£dx(ad-L,.,X,)(£d,(L,„ 

= [[(Af)*X2,X,]]. 
The other recursion operators are treated analogously. Thus we obtain: 
(Af)* = Af, (Af)*=A^. 



6 Fundamental properties of the NLEE 
6.1 Integrals of motion 

In this subsection we are going to apply a method by Drinfel'd and Sokolov [7] to derive the 
integrals of motion to system (2.1). In order to do that it proves to be technically more convenient 
to deal with the Lax pair (2.13), (2.14). Then we map the operators L and A into 

C = V-^LV = id^ + \J + Cq + ^ + --- , (6.1) 

A 

A = V-^AV = idt + + \A-i + A + 4^ + • • • , 

A 

where all matrix coefficients £fc, A-i and Ak, k = 0, 1, . . . are diagonal and make use of the 
following asymptotic expansion for V{x,t, A): 

P(x,t,A) = IL + ^ + ^ + .... 

A A^ 

In order to make all further considerations unambiguous we assume that all coefficients pi 
(Z = 1,2,.. .) are off-diagonal matrices. 

The zero curvature representation is gauge invariant, i.e. [C,A] = is fulfilled. Since 
= the commutativity of C and A is equivalent to the following requirements 

dxA^i = 0, dtCk - dxAk = 0, A; = 0, 1, ... . 

Hence Ck represent densities of the integrals of motion we are interested in. 
Equality (6.1) rewritten as 

Lv = vc 

holds identically with respect to A. This requirement leads to the following set of recurrence 
relations: 

Uo + Jpi = Co+piJ, 
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Figure 6. Plots of the density of the first integral of motion as a function of x evaluated on the 
stationary quadruplet soliton (3.26) ioia — I3 = j = uj^l, 6^0 (left panel) and for the stationary 
doublet soliton (3.28) for k = do = 1, = (right panel). 

ipi,x + Uopi + Jp2 = Ci+ piCo + P2J, (6.2) 

k-l 

ipk,x + Uopk + Jpk+i = Ck+Pk+iJ + ^ 

m=0 



In order to solve it we apply the well-known procedure of splitting each relation into a diagonal 
and off-diagonal part. For example, treating this way the first relation above one gets 



f 



[J, Pi] 



(6.3) 



'0' •-'0 

where the superscripts d and / above denote projection onto diagonal and off-diagonal part of 
a matrix respectively. Taking into account the explicit form of Uq (see formula (2.15)) for Cq 
we have 

10 



£0 = -{uul + vvD 



0-2 
1 



Thus as a density of our first integral we can choose: Iq = u*Ux + v*Vx- It represents momentum 
density of our system. For the stationary solutions (3.26) and (3.28) the momentum density is 
depicted on Fig. 6. It is evidential that for both cases of solutions it is a well localised function 
of X. On the other hand after inverting the commutator in the second equation in (6.3) one 
obtains 

V^{uvx-vux) {uu*+vv*)/2 

p, = - adf ui = -'^\ V2{u*v*x - v*u*x) -V2iu*v*x - v*u%) \ .(6.4) 



{u*Ux + v*Vx)/2 



"X 





Similarly, for the second matrix Ci one needs to extract the diagonal part of (6.2). The result 
reads 



{uIp,Y = c,. 

After substituting the expression (6.4) for pi into (6.5) one obtains 

£1 = - + VV*x\^ + 4:\uVx — VUx\^\ 

^ V ( 

Hence the second integral density is Xi = \uu*x + vv*x\'^ + ^\uvx — vux\^- 



(6.5) 
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In general, one is able to calculate the-Zc integral of motion through the formula 

The matrix pk in its turn is obtained from the following recursive formula 
/ fc-i \ 



Pk 



a,d/ ipk-i,x + {UoPk-iY - ^ Pk—l—m 



m=0 



6.2 The hierarchy of NLEE — Lax pair approach 

Our goal here is to describe the hierarchy of integrable NLEEs associated to the Lax opera- 
tor (2.3) in terms of the operators Ai and A2. In order to achieve it we shall analyze the 
recurrence relations obtained form the zero curvature representation [L, A] = 0. 

An arbitrary member of the integrable hierarchy under consideration has a Lax pair in the 
form 

N 



L{\)=id.^ + XLiix,t), A{X)=idt + J2^''M^,t)- 

k=l 

As before the operators L and A are subject to the action of the reductions (2.7), (2.8). Hence 
the coefficients of A are hermitian matrices to fulfill 

CA2,_lC = -^29-1 ^ ^2,-1 GS[1(3), CA2gC = A2g ^2gGs[0(3). 

Thus the original Lax pair (2.3), (2.4) represents the simplest nontrivial flow {N = 2) of the 
above general flow pair. 

The compatibility condition [L(X), A{X)] = gives the following set of recurrence relations: 

[Li,An]=0, (6.6) 

id^AN + [Li,AN-i]=0, (6.7) 



id^Ak + [Li,Ak-i] = 0, k = 2,...,N-l, (6.8) 



d,Ai - dtLi = 0. (6.9) 

It directly follows from the first relation that the highest order term is a polynomial of Li. Since 
Li G s[^(3) and L2 G s['^(3) we have two options for A^- 

a) AN = f2pL2 for N = 2p, b) An = f2p+iLi for N = 2p+1. 

It suffices to restrict ourselves with the case when N = 2p since the case N = 2p+l is completely 
analogous. 

As we did many times in our exposition we shall split each element Af^. into two mutually 
orthogonal parts A2q_i and f2q-iLi (resp. A-^q and f2qL2 in the case even indices): 

A2q-l = A^.i + /2g-l^l, = + /2q^2. (6.10) 

Substituting the splitting of A^-i into (6.7) we have 
if2p,xL2 + if2pL2,x + [Li,A2p_i] = 0. 
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After taking the Killing form ( • , L2) to separate the Li-commuting part and its orthogonal 
complement we deduce that fN = CN = const and 

^2p-i = -^C2p ad^^^ L2,x- 
Similarly, after inserting the splitting (6.10) into a generic recurrence relation (6.8) we obtain 

if2q-l,xLl +if2q-lLl,x +^(^2g-l)^ + [^l>-42g_2] = 0, 
if2q,xL2 + if2qL2,x + ^ (^^q) ^ + = 0. 

After extracting the Li-commuting part from the equations above one obtains the coefficient 

f2q-l (resp. f2q) 

1 3 

f2q-l = C2q-1 " ^ ^ (^^g-l) ^ ' ^l) . = C2q - l^d-\{A^g) ^, L2) , 

where C2q-i (resp. C2q) is a constant of integration. On the other hand for A2g_i (resp. A^^) we 
have: 

= ^1 (^2g+l) " «C2g+l ad^^^ Li^x, -4^g-l = A2(^2g) " ^C2g ad^^^ L2,x, 

where the integro-differential operators A2 and Ai are given by (4.10) and (4.11) respectively. 
The last recurrence relation (6.9) yields to 

/i = ci - ^d^\{Ai)x,Li), iadll dtLi + AiAi - ici ad^l Li,x = 0. 

What remains is to substitute consequently the expressions for Aj:, k = 1, . . . , 2p — 1 above in 
order to write the NLEE in terms of the operators Ai and A2. Here we give the results for both 
cases a) and b): 

p p-i 

a) ad;^;^^ dtLi - C2g(AiA2)''"Ui ad;^^^ L2,x - J] C2g+i(AiA2)« ad;^^^ Li,^ = 0, 

P P 

b) ad^^^ dtLi - C2g(AiA2)'?-Ui ad;^^^ L2,x - ^ C25+i(AiA2)'? ad^^^ Li,^ = 0. (6.11) 

g=l q=0 

The coefficients are involved in dispersion laws of NLEEs. By analogy with (2.18) the 
dispersion law of a general NLEE is defined through 

N 

/(A)= hm g^'T^^M^,t)9±- 

x^±oo ' 

k=\ 

The dispersion law governs the time evolution of scattering matrix (2.19) through the linear 
equation 

iatr(t,A) + [/(A),r(t,A)] = ^ r(t, a) = e*-^(^)*r(o, A)e-*-^(^)*. (6.12) 

It is not hard to check that the equalities below are valid 

a) /(A) = Y C25+iA'''+ V + Y ^2'?^"'^' 

g=0 q=\ 

b) /(A) = C2.+iA'^+V + Y ^25^''^- (6-13) 

g=0 q=\ 

We remind that the constant diagonal matrices J and / represent the diagonal forms of the 
asymptotic of L\ and L2 respectively, see (2.13), (2.14). As one can easily convince himself the 
initial NLEE (2.1) can be derived from the above formulae in the simplest case N = 2 after 
plugging C2 = — 1 and ci = 0. 
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6.3 The hierarchy of NLEE and generahzed Fourier transforms 

In the previous subsection we described the class of NLEE related to the Lax operator L. In 
particular we showed that if Li{x,t) is a solution to one of the NLEE (6.11) then the corre- 
sponding scattering matrix T{X,t) satisfies the linear evolution equation (6.12) with dispersion 
law /(A) given by (6.13). 

Here we will prove a more general theorem. 

Theorem 1. Let the Lax operator L be such that its potential satisfies the conditions Cl)-C3). 
Then the NLEE (6.11) are equivalent to each of the following set of linear evolution equations: 

.dS+ 



dt 



+ [/(A),S+(A)] = 0, 


dD+ 
' dt 


= 0, 


.dS- 

' dt 


+ [/(A),5-(A)]=0, 


+ [/(A^)>4]=0, 


dt 


= 0, 


dt 


+ [/(Afc),^,-]=0; 


+ [/(A),r-(A)] = o, 


dD- 
' dt 


= 0, 


dT+ 
' dt 


. + [/(A),T+(A)] = 0, 


+ [fiXt),r^] = 0, 


'dt 


= 0, 


M 

dt 





dt 

^^ + [/(At),4]=0, ^^ = 0, i"-^ + [f{X~),s-,]=Q; (6.14) 

dT- 



Proof. Let us first consider the NLEE (6.11a) and let us denote its left hand side by Ng^-^{x,t). 
Next let us expand it over the complete set of eigenfunctions H^^{x,X). If we make use of 
the expansion (5.11) then we need to calculate the expansion coefficients (5.12). Using the 
properties of the recursion operators (5.16) and (5.17) we have: 

7l)(A)= [[e^a(y,A),ad^iiV^)(y)]] 

= [[i/±Jy,A),Li,i]] -f]c2,[[if|„(y,A),(A-'*)'?-iArad^iL2,.]] 

9=1 

p-1 

- J2 C2g+i [[kUv, A), (A-'*)^ ad^i L,4 

9=0 

= ^<5ir±(A) + f]c2,[[A+(A+)^-ii^±Jy,A),ad-iL2,.]] 

9=1 

p-1 

+ J2 C2,+i [[{A'-yK^Ay, A), ad^i Li,,]] 

9=0 

= {^tr^W + E ^29^'^-' [[H$M, A), ad^i L2,.]] 
9=1 

p-1 

+ C2,+iX^'+' [[K^M, A), ad^^i Li,,]] 

q=0 

= j(^i(^S^^,E^a'^ + {S^f{X)S^{X),E^c.?j=0 for all a G A+. (6.16) 
From equation (6.16) we find that 

f)Q± 

^S^— + S+f{X)S^{X) = n{X), (6.17) 
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where 'H(A) G [}. It remains to take the Killing form of (6.17) with the Cartan elements J and / 
of si to determine that ^(A) = /(A). Thus we have proved that from the NLEE one gets the 
first of the equations in (6.14): 

f)Q± 

Next we rewrite equation (2.23) in the form: 
D+f-T+D-{\) = S+S-{\), 
and equate the principle upper- and lower-minors of both sides. As a result we get: 

1 _|_ _j_ 1 . _|_ _j_ _j_ , 1 _|_ _j_ 

+ - ~ 1 ~ ^ai^a-i ~ '^03 '^03 + ^(•^ai'^a2'^a3 ~ '^ai '*«2 '^03 ) '*a2 '^"2 ' 

TTl-^ T7I2 ^ ^ 

+ — ^a2^a2 ~ ^a-i^a-j, ~ 9 (-^ai '^02 '^"3 ~ '^ai'*02^a3) '^oi ^^2 '^"2 (6.18) 

and 

1 !Y* f (y* y>» I fy* ( ty I ty^ I fv* fv* ly^ rv* I ^ I ry* I ly* ly^ I ly* 

+ — ^ ' a2 a2 ' as' as ai' a2 as ' ai' a2 as' ~^ a ai' ai' a2 a2'' 

1 /y* ^ /y> ry» I /va I f fy* I /y^ I ty ty ry ly I ^ I ry I ry ry I ry { f~\ 1 

-(- — ^ ' cii ' Qi ' a^' a-j ^ r\\' ai' a2 ct'i ' ai' 0.2 I ' a.\' ai' 0.2 02 yj.i^^ ) 

for A G M. We remind that the scalar functions mf{X) and m^(A) (resp. m^(A) and rrig (A)) 
are analytic functions for A € C+ (resp. A G C_). Then the relations (6.18) can be viewed as 
a RHP which can be solved by the Plemelj-Sokhotsky formulae. Thus equations (6.18) allow 
us to recover mf{X) and m^(A), and effectively D^{X) in their whole regions of analyticity 
from S^{X). Likewise, equations (6.19) allow us to recover D^{X) from T^(A). 
From equations (6.18) there also follows that 

—A- = k = l,2, i.e. —— = 0. 
dt dt 

We can also reconstruct T^(A) as the Gauss factors of D~ (X) and check that 

idtT^it,X) + [fiX),T^{t,X)] = 0. 

Thus we have proved the statement of the theorem for the scattering data on the continuous 
spectrum of L. Similar procedure allows us to recover the equations for the data on the discrete 
spectrum of L. 

In order to complete the proof of the equivalence between the NLEE (6.11) and (6.14) it 
remains to apply Corollary 1. 

The equivalence of the NLEE (6.11) and (6.15) is proved along the same lines analyzing the 
expansion coefficients 6^ ^(A). ■ 

As an immediate consequence of the theorem we obtain: 
Corollary 2. Each of the minimal sets of scattering data S^, k = 1,2: 

51 = {4(A),s-(A), AgM, <fc,s;.fc,A±, aGA+, = 1, . . . , 2iVi + A2}, 

52 = {r^(A),r-(A), A G M, r+.„r-„A±, a G A+, k = 1, . . . ,2Ni + N2} , 

determines uniquely both the scattering matrix T(X) and the potential Li{x) of the Lax opera- 
tor L. 
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Proof. The fact that each of the sets determine the scattering matrix T(A) foUows easily 
from the Theorem. Indeed, we showed how, starting from Sk one can construct each of the 
Gauss factors of r(A); so it remains just to take the corresponding products. 

In order to reconstruct the corresponding potential we have to solve the RHP with canonical 
normalization for A = 0. Considering the Taylor expansion of A): 

oo 

x^(x,A) = a + j;A^xf,)(x), 

n=l 

we get 

Li[x) = -I — y — . ■ 

ax 

6.4 Hierarchy of Hamiltonian formulations 

Our analysis here is again based on the Wronskian relations, extending the ones in Section 4. 
We first use slightly modified equation (4.18): 



{x^6x^{x,X),J)\'^^ = iX / dx{6L,,h^{x,X)) 

J —oo 

/oo 
dx {6L,,hfix,X)), 



where 

h^{x, A) = x^{x, X)Jx^{x, A), hf{x, A) = x^{x, A)/x^(x, A). 

A third class of Wronskian relations connects the variation of the derivative xi^, A) 
dxx{x,X) with Li, i.e. 



{x^X^{x,X),j)\'^^ = - / dx{Li,h^{x,X)), 

J — oo 

/oo 
dx {L,,hfix,X)). 



Thus we derive the interrelations between the generating functionals of integrals of motion 
and the corresponding potential: 

{D^6D^,J) = ±6ln{mfmf) = -iX[[adll 6Li, h'^ix, X)]], 

{D^6D^,I) = 6ln'^ = -iX[[adll6Li,hf{x,X)]]. 



Note that the second reduction in equation (2.21) applied to /ij(x,A) and hf{x,X) gives 
Ch^{x, -A)C = -h^{x, A), Chf{x, -A)C = hf{x, A), 
which means that 

oo + oo 

mi 



where Cj-2s+i and Ci-2s are integrals of motion. 
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Let us note that neither h^{x, A) nor hf{x, A) are eigenfunctions of the recursion operators. 
Indeed, if we separate the 'orthogonal' and 'parallel' to Li and L2 parts and split them into: 



/t±'^(x,A) = Hj{x,X) + Kj{x,X), hf'^{x,X) = Hf{x,X)+Kf{x,X), 



we get: 



AfKj{x,X) = XHj{x,X) + iadllLi^x, AfHj{x,X) = XKj{x,X), 

AfKf{x,X) = XHf{x,X), AfHf{x,X) = XKf{x,X) + iadllL2,x. 

Therefore 

(Af Af - X^)Kj{x, A) = iAf ad^i Li,,, (Af A± - X^)Hj{x, A) = iXadl^ Li,,., 

(Af Af - X^)Kf{x, A) = fAad;^^^ L2,x, (Af - X^)Hf{x, X) = iAf adl] Li,^ 

and 

00 

±51n(mf mf ) = ^ X-^'-HCj,2s+i = A[[ad^i JLi, (A2A1 - X^)-^A2 ad^^ Li,,]] , 



s=0 

-|- 00 

<51n = ^ X'^'5Ci,2s = X^ [[ ad^i 5Li, (A2A1 - A^)-! ad^^ , 

^ S = 



or in other words 

5Cj,2s+i = - [[ad^i (A2Ai)^A2 ad^^ (5Li,,]] , 

'^C/,2. = - [[ adl\ 5Li, (A2Ai)^ ad^i 5L2,4 . (6.20) 
We can rewrite equations (6.20) in the form: 

Thus in fact we have proved that the conserved quantities of these NLEE satisfy the well known 
Lenart relation 

oLi[x) dLi[x) dLi[x) oLi(x) 

We finish by formulating the Hamiltonian properties of these NLEE. It is only natural that 
the Hamiltonians must be linear combinations of integrals of motion: 

p-i p p p 

Ha) = ^ C2q+lCj^2q+l + ^ C2qCl^2q, -f^b) = X] C2q+lCj^2q+l + ^ C2qCl^2q- 

q=0 q=0 q=0 q=0 

Next we introduce a symplectic structure using the symplectic form: 
= [[ ad;^;^^ 6L1 A A2 ad^^^ 5Li\] . 

Using (6.21) we find that the equation of motion (6.11) can be written down as 

Q\ = [[ adll 6L1 • A2 ad;^^^ = 6H, (6.22) 
with H = iJg,) or H = i^t,) respectively. 
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In particular, choosing H = Cj^2 we find that equation (6.22) becomes 

and coincides with the reduced HF equation (2.1) we started with. 

Along with Qi we can also introduce the following hierarchy of symplectic forms: 

n2k~i = [[ ad^i 6Li A(A2Ai)^-U2 ad^i 5Li]] . (6.23) 

Then we can write the equation of motion (6.11) in Hamiltonian form using each one of the 
above 2-forms as follows: 

f^2fc-i= [[ad^i^^i(A2Ai)'-'A2ad^iLi,i]] = 6H2k 

with H = i^2fc;a) or H = i/2fc;b) respectively, where 

H2k;a.) = C2g+2fc+l<^J,2g+l + C2qC I ^2q+2k , 
q=Q q=0 
P P 

-f^2fc;b) = C2q+lCj^2q+2k+l + C2qC J ^2q+2k ■ 

q=0 (j=0 

Using the expansions of ad^^^ 5Li over the 'squared solutions' 

^2k-i = [[ad-ijLi A(A2Ai)*^-iA2ad-i5Li]] = — / d\\^^-\^+ {X) - {\)) 

2N1+N2 / \ 
+ y Res \^^-^n^{X) + Res X^'^^QiX) , 

where 

^tW= E '^V^(A)A5'r±(A). 
Thus we conclude that the family of symplectic forms (6.23) are dynamically compatible. 

7 Discussion and conclusions 

A system of coupled equations, which generalize Heisenberg ferromagnet equations has been 
studied. The system is associated with a polynomial bundle Lax operator L related to the 
symmetric space SU{3)/S{U{1) x U{2)). 

The spectral properties of the operator L in the case of the simplest constant boundary 
condition (2.2) have been described. The continuous spectrum of L fills up the real axis in 
the complex A-plane and divides it into two regions: the upper half plane C+ and the lower 
half plane C_. Each region is an analyticity domain of a fundamental analytic solution to the 
auxiliary linear problem. 

By using the dressing method [40, 49, 50, 51] we have constructed 1-soliton solutions of the 
(Z2-HF) models. Due to the additional symmetries the Lax operator L may possess two types 
of discrete eigenvalues: generic ones and purely imaginary ones. Therefore we will have two 
types of soliton solutions - quadruplet and doublet solitons respectively. We outlined the purely 
algebraic construction for deriving the A/'-soliton solutions for both types of solitons. 



Polynomial Bundles and GFT for Integrable Equations on A. Ill-type Symmetric Spaces 43 



Using the Wronskian relations one is able to construct 'squared solutions' and an integro- 
differential operator called recursion operator whose eigenfunctions they are. There exists an- 
other viewpoint on recursion operator - they generate hierarchy of symmetries of NLEEs. Thus 
one can derive the recursion operator of a NLEE from purely symmetry considerations. 

Using the interpretation of the ISM as a generalized Fourier transform and the expansions over 
the 'squared solutions' we studied the fundamental properties of the class of NLEE admitting 
a Lax operator of the form (2.3) , including the description of the whole class of NLEE, the infinite 
set of integrals of motion and the Hamiltonian formulation of the corresponding hierarchy. 

The results of the present article can be extended in several directions. Firstly, one can 
develop the theory in the case of a rational bundle L: 

L = idx + ALi + 

A 

The simplest NLEEs related to Lax pair of this type take the form 

iUt = Uxx - {u{u*Ux + V*Vx))x + 8vv*u, 

iVt = Vxx - {v{u*Ux + V*Vx))x - 8uu*v, (7.1) 

where u and v are functions of x and t subject to the same condition as in the polynomial case, 
i.e. |up + \v\'^ = 1. The system (7.1) can also be seen as an anisotropic deformation of (2.1) with 
k = 1, N = 3 and is a special case of the models proposed in [23]. 

That modification is required when one imposes an additional Z2 reduction of the form 
A — )• 1/A, see [18]. This case is more complicated and much richer than the one we have studied. 
It requires the construction of automorphic Lie algebras and studying their properties following 
the ideas of [30, 29, 46]. 

The second direction of generalization concerns considering Lax operator L related to a ge- 
neric symmetric space of the type A. Ill = SU{n + k)/S(JJ{n) x U{k)) or, more generally, 
related to other types of symmetric spaces. This will allow one to treat various multi-component 
generalizations of our NLEEs. Geometric properties of similar models related to other symmetric 
spaces are studied in [47]. 



A Parametrization of Gauss factors 

Here we will list useful formulae allowing one, given the scattering matrix T(A) to evaluate its 
Gauss factors 5^ (A), T^(A) or rather s^(A), r^(A) 

S^{\)= exp (A) , (A) = exp (A) , 

st, st,\ / r+ r+ 

s+(A) = I 4 1 , r+(A) = r+ 

000/ \000 

ooo\ /ooo 

s-{\)= \ s-_^ , r-(A)= r^^ 

D^{X) = d\a.g{ml ,1712 /m'l , l/"^2^)i ^ W — diag(l/m2 , m-j m^). 

Next we insert these formulae in the left hand sides of the Wronskian relations (4.21), (4.22) 
and (4.23). Skipping the details we find 

-'12 + -'l3-'21 — ^23^11 + -tl2-'23 — -t22-tl3 -^13 
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Table 1. The expansion coefficients of ad^^^ ij.xi ^'^lI ^-^i ^"^^ ^'^lI ^i.t- 



Roots 


a\ = e\ — 62 


ai = 61 — 62 


ai = 61 — 62 


(a, J) 
(a, I) 


1 
1 


1 

-1 


2 





<(A) 
-V(A) 


4, (A) 

-■Sa2(A) 


243(A) 
-25^3 (A) 


rP--(A) 


<(A) 

-Sai(A) 


-<(A) 

Sa2(A) 


4i 42(A) 

^ai ■5a2(A) 




(™^)^^+ 


("^2^)^^+ 

(A) 


-2m^m^r+g(A) 




H^Vi(A) 

(^2)^- (A) 


^""^^'s- (A) 


m5^m+r+ s+ (A) 

777,^1 777,1 r~ ,S~ I A 1 




5<(A) 
.5s- (A) 


542(A) 
5Sa2(A) 


5s+(A) + i5<.+ (A)-i5s+<(A) 
5sa3(A) - ^55-5+ (A) + ^^s-^s" (A) 


S'tr-{X) 


^4i(A) 
|V(A) 


l<(A) 

S'5a2(A) 


3t43(A) + 2'^a2(A)5i^ai ~ 24i(A)3j42 
S'^a3(A) ~ 2^a2(A)at^ai 2 '^ai (A) S'^a2 



^21-132 — J22J31 J31 

IE IlT 5 

-t 33-^21 — ^22^31 _ -^31 

13^22 — 23-^11 _|_ ^13 
2m^ 2m^ ' 

{S+JS+{X),E^) = (a,J)4(A), 

(S'+/S'+(A),£'q,) = (Q,/)s;i;(A) for = 61-62,62-63, 

(5+/5+(A),i?ei-e3) = 4243; (A-l) 

(5+55+ (A), E'a) = (54(A) for = 61-62,62-63, 
(5+55+(A),i?e,-e3) = SS+ + ^(5(4243). 

Thus we are able to express all expansions coefficients in terms of s^. The results are: 
r«'±(A) ^ (5±J5±(A),i?^„) = ±(a, J)4(A), 

42).±(A) = {S^IS^{X),E^^) = {(±(a, J)4(A) for a G (5+)} ± (a, J2)4(A). 

Note that varying Li may lead to varying the positions of the discrete eigenvalues. In order 
to evaluate the corresponding expansion coefficients, now we make use of the Wronskian rela- 
tions (4.23). We use also the condition C2 and the explicit formulae that allow one to express 
all the Gauss factors 5^ (A), T^{X) and D^{X) as rational expressions of the matrix elements 
of r(A). Condition C2 ensures that the Gauss factors have at most first order poles in the 



-t21 _ -t 11 -t 32 — J 12 J 31 

ai + ' a2 + ' as 

- _ -L32 - _ -t 23-^31 ~ ^21-^33 

^a2 —1 ^ai — ' ^03 

1712 



-t23 4. -'l2-'33 — -tl3J32 

h 



' a2 — ' ai — 5 '03 
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vicinity of A^. Therefore the variations of S^{X) and T^(A) in the vicinity of take the form: 
Ss^iX) (A - A^)^-+i {-(X^^^^^ + + ^^'^) ' 

where s^.^ and r^.^ determine the values of the corresponding Gauss factors for A = A^. Thus, 
along with the (4.23) we obtain: 

Pik^^i = [[eL;„ adZ^i 6L,]] , 6' p%, = i [[e%^.^, + A±e|,^„ ad^^^ 5L,]] , 

Remark 3. More precise treatment of the contribution from the discrete spectrum starts by 
considering the potentials for which Li[x) — Li_as is on finite support. Then the Jost solutions 
of -L, as well as the scattering matrix T(A) become meromorphic functions of A which ensures 
the validity of all considerations above. The next step would be taking the limit to potentials 
of Schwartz- type. These considerations come out of the scope of the present paper. 

When we analyze the discrete spectrum we make use of the conditions C2 and C3. Since 
the discrete eigenvalues are zeroes of the principal minors, and assuming that A^ (resp. A^) is 
a zero of, say m'^{X), then we have: 

mf(A) = (A-A±)m±, + 0((A-A±)2), 

and similar expansions for m^(A). Using remark 2 on formulae (A.l) we easily find, that 
Ta^'^{X) have simple zeroes in the vicinity of A^. 

We end by reminding a few useful formulae of frequently appearing expressions in the main 
text; they are derived in [18, 19]. 

Since Li satisfies the characteristic equation L\ = Li then ad^^ has eigenvalues ±2, ±1, 
and satisfies the characteristic equation: 

(adi^ -4)(adi^-l) adi, =0 
and therefore 

ad^i = ^(Sadii -adij. 
If we choose A S sl^ we have 



ad^;^X=( u*a-a*u u*b-va* \ + -{w - w*) 

v*a — ub* v*b — b*v J 

where w = ua* + vb* . If in particular, X = Li^x we find: 

















U V 


v*u 





l\wQ 

£1 = adj^_^^ Li^x = u*Ux — u^u + |i(;o|iip u*Vx — wn* + ^■wou*v 
\ v*Ux - uvl + ^wquv* v*Vx - v*v + Iwolvl"^ 

where wq = uu^ + vv*. Similarly for Y G sf and ad^^^ Y we have 

-k — n \ ^ / uao + vm* vai + um 

Y = \ k m \ , adll ^ = ^ -u*ao - mv* 

m* n / \ —v*a^ — u*m* 
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aQ = 5k + n- 3W, ai = k + 5n - 3W. 
Choosing Y = L2^x we get: 

(0 —u\u\'^ — v{uv*)x — — u(f 

u*\u\l + v*{u*v)x 

v\v\l + u{u*v)x 

—Ux — u{uu*^ + vv*) —Vx — f (nn* + vv*) 

ul + U*{u*Ux + V*Vx) 

vl + v*{u*Ux + v*Vx) 
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